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FOREWORD 


This  report  examines  several  analytical  and  instrumentation  methods 
for  eliminating  the  interferences  caused  by  reflection  of  sound  waves 
from  the  ground  surface  during  FAR  Part  36  aircraft  noise  certification 
measurements. 

The  author  acknowledges  the  significant  contributions  made  by  C. 
Bartel  for  investigating  the  aircraft  position  accuracy  problem,  D.  Hoy 
for  researching  surface  impedance  data  and  programing  several  key 
computer  modules,  R.  Helizon  and  A.  Segal  for  aircraft  noise  data 
acquisition  at  Los  Angeles  International  Airport,  L.  Sutherland  for 
technical  guidance,  and  J.  Parkinson  for  overall  program  management. 

This  work  constitutes  Technical  Areas  1  and  5  under  Contract  DOT- 
FA7  8WA-4143 . 
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0.  INTRODUCTION 


Administration  of  the  rules  embodied  in  FAR  Part  36  -  Noise  Standards: 
Aircraft  Type  and  Airworthiness  Certification,  has  involved:  (j)  adherence  to 
specific  procedures  contained  in  Appendices  to  the  regulation,  and  (2)  application 
of  administrative  processes  to  approve  or  define  those  test  procedures  not 
specifically  called  out  in  the  regulation.  It  is  this  latter  aspect  of  administering 
FAR  Part  36  that  this  report  addresses.  FAA  staff  in  Washington,  and  in  regional 
offices,  have  developed  many  general  and  specific  guidelines  not  codified  within 
existing  Federal  Aviation  Regulations.  However,  the  guidelines  for  noise  certifi¬ 
cation,  such  as  those  provided  in  the  handbook  on  certification  procedures 
published  by  the  FAA  Western  Regional  Office*  do  not  include  explicit  details  on  a 
number  of  procedures  such  as  the  following: 

o  Methods  for  correcting  measured  noise  spectra  for  ground  reflection 
effects 

o  Methods  for  correcting  measured  levels  when  ambient  noise  levels  are 
excessive 

o  Methods  for  correcting  propeller  aircraft  noise  data  for  deviations  in 
aircraft  performance  from  reference  conditions 

While  informal  rules  have  generally  been,  established  by  the  FAA  out  of 
necessity  for  each  of  these  correction  procedures,  they  have  not  yet  been  fully 
documented  or  validated. 

The  objective  of  this  document,  Volume  I  of  a  multiple  volume  set  on 
correction  procedures,  is  to  provide  the  technical  background  required  to  allow 
FAA  to  codify  ground  reflection  correction  procedures  in  future  revisions  to  either 
administrative  guidelines  handbooks  or  to  FAR  Part  36  itself.  In  the  latter  case, 
statements  in  the  existing  regulations  allowing  "approved  procedures"  or  "approved 
corrections"  could  be  replaced  with  explicit  details. 

FAR  Part  36  data  correction  procedures  can  be  categorized  into  two  general 
classes: 
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o  Procedures  that  correct  for  nonstandard  test  conditions  (e.g.,  correc¬ 
tion  for  sound  attenuation  at  nonstandard  temperature  and  humidity)  or 
aircraft  performance. 

o  Procedures  that  compensate  for  limitations  of  the  measurement  and 
sound  description  methods  (e.g.,  correction  for  aircraft  sound  levels 
masked  by  ambient  noise  or  reflections  off  the  ground). 

Corrections  for  nonstandard  performance  of  propeller  aircraft  fall  in  the  first 
category  and  corrections  for  ground  reflection  or  ambient  noise  levels  fall  in  the 
second  class.  Table  0-1  identifies  the  major  types  of  data  corrections  that  are 
explicitly  specified  in  FAR  Part  36.  Table  0-2  identifies  these  data  and  procedures 
which  are  not  completely  specified  and  which  require  FAA  approval. 

Applicable  paragraphs  identified  as  "A,"  "B,"  or  "C"  in  Tables  0-1  and  0-2 
refer  to  large  subsonic  transport  airplanes  and  all  turbojet  powered  aircraft 
regardless  of  category.  Paragraphs  identified  as  "F"  refer  to  light  propeller  driven 
aircraft  (A,  B,  C,  and  F  are  Appendices  within  FAR  Part  36).  Paragraphs  identified 
with  an  (*)  are  associated  with  the  technical  correction  procedures  of  this  series  of 
volumes.  This  volume  I  addresses  only  the  pseudotone  (ground  reflection)  problem. 
At  the  time  of  writing,  other  volumes  concerning  ambient  noise,  propeller  aircraft 
performance,  and  appropriateness  of  tone  correction  were  either  in  preparation  or 
being  contemplated. 

The  significant  features  of  this  volume  are  as  follows: 

o  Analytical  methods  of  pseudotone  removal  are  not  practical  at  this 
time. 

o  The  10  m  (pole)  microphone  position  is  recommended  for  certification 
measurements.  This  is  expected  to  effectively  eliminate  pseudotones 
at  takeoff  and  approach  locations  without  any  further  correction. 
However,  the  800  Hz  low  frequency  cutoff  for  tone  corrections  should 
be  retained  for  the  sideline  position  only. 

o  A  substantial  number  of  aircraft  flyover  noise  measurements  were 
obtained  in  support  of  this  study.  These  are  summarized  in  this  volume. 


Table  0-2 


Data  or  Procedures  Requiring  FAA  Approval 
Currently  Specified  in  FAR  Part  36  (Change  7) 


Category 

Data  or  Procedure 

Applicable  Paragraph 

Aircraft  Performance 

Aircraft  Position  Measurement 

A36. 1(d) 

A36. 1 1 

F36. 109(f)  (6) 

Sample  Rate  for  Aircraft 
Performance  Data 

A36. 1(d)  (1) 

A36.3(b) 

F36. 109(g) 

Aircraft  Reference  Flight 
Profile,  Speed,  or  Power 

Setting 

A36.ll 

C36.7(f) 

C36. 9(f)(1) 

F36.ll 1(b)  (1) 

Aircraft  Performance  Measure¬ 
ment  Instrumentation 

A36.5 

F36. 101(b)  (6) 

Aircraft  Noise 

Corrections  for  Non-standard 
Weight 

a 

A36. 1  (7) 

Corrections  for  Non-standard 
Flight  Profile  Speed  or  Power 
Setting 

A36.ll 

Noise  Measurement 

Equipment 

A36.3 

F36.103 

F36.105 

Ambient  Noise  Corrections 

A36.5 

F36. 107(c) 

Overall  Data  Corrections 
Procedures 

A36.3 

A36.5 

F36. 109(a) 

F36.203 

Atmospheric 

Para  me  ten 

Measurement  Facility 

A36.3 

A36.5 
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I. 


GROUND  REFLECTION  AND  PSEUDOTONES 


l.l  Nature  of  Ground  Reflection  Phenomenon 

Technical  Area  I  of  this  document  is  concerned  with  the  influence  of 
sound  reflection  from  the  ground  surface  on  measured  aircraft  noise  spectra  and 
resulting  EPNL  values.  The  word  "pseudotone"  has  become  common  to  describe  the 
spectral  irregularities  caused  by  the  ground  reflections  in  the  lower  frequency 
portion  of  the  spectrum  (generally  below  I  kHz). 

The  measurement  setup  required  by  FAR  Part  36  is  sketched  in  Figure 
l-l.  In  addition  to  the  direct  wave  from  the  aircraft  to  the  microphone,  the 
indirect  wave  reflected  from  the  ground  surface  is  also  sensed  by  the  microphone. 
Direct  and  reflected  waves  instantaneously  add  which  leads  to  an  interference 
pattern,  i.e.,  a  succession  of  constructive  and  destructive  interferences.  The 
nature  of  this  pattern  depends  on  the  geometrical  relations,  the  ground  surface, 
and  the  sound  frequencies.  Figure  1-2  shows  an  example  of  an  interference  pattern 
in  the  form  of  a  theoretical  correction  chart  for  the  case  of  source  and  receiver  at 
the  same  height  over  an  infinite  impedance  ground  surface.  The  dB  correction  is 
subtracted  from  a  measured  spectrum  in  order  to  obtain  the  spectrum  that  would 
have  been  measured  under  free  field  conditions.  Although  this  is  not  typical  for  a 
flyover  application,  the  figure  shows  several  features  of  general  interest.  At  the 
lower  frequency  end,  the  correction  approaches  6  dB  which  corresponds  to  the 
well-known  doubling  of  acoustic  pressure  near  a  perfectly  reflecting  surface.  At 
the  upper  frequency  end,  the  correction  for  discrete  tones  consists  of  very  many 
closely  spaced  lobes.  This  is  not  representative,  however,  since  spectrum  analysis 
of  aircraft  noise  is  performed  with  finite. bandwidth  filters.  The  wider  the  filter 
bandwidth,  the  faster  the  correction  will  tend  towards  the  stable  limit  of  3  dB 
corresponding  simply  to  the  energy  summation  of  two  uncorrelated  signals  of  equal 
strength.  For  the  discrete  frequency  curve,  destructive  interference  leads  to  an 
infinitely  large  dip  for  a  pure  tone.  For  a  finite  bandwidth  of  noise,  such  an 
infinitely  large  dip  cannot  occur;  the  wider  the  bandwidth,  the  less  pronounced  this 
dip. 

An  example  of  pseudotones  observed  during  a  real  flyover  is  shown  in 
2 

Figure  1-3.  A  series  of  spectra  at  regular  time  intervals  are  shown.  The  theory 
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Reflection  Correction:  Point  Source,  Infinite  Impedance  Boundary 


of  ground  reflection^  shows  that  interference  pattern  frequencies  are  inversely 
proportional  to  the  path  length  difference  between  the  direct  and  reflected  paths. 
It  is  easy  to  see  that  this  difference  is  largest  when  the  aircraft  is  overhead,  and 
decreases  as  the  aircraft  is  at  smaller  elevation  angles.  The  pattern  in  the  figure 
therefore  explains  itself:  pseudotones  at  higher  frequencies  during  the  approach, 
decreasing  to  lower  frequencies  around  over  head,  and  increasing  again  as  the 
aircraft  recedes.  This  is  in  marked  contrast  to  the  behavior  of  the  Doppler  shift  of 
true  (fan)  tone  components  also  indicated  in  the  figure:  it  starts  at  an  elevated 
frequency,  experiences  the  most  rapid  frequency  change  around  overhead,  and 
settles  at  a  lower  frequency.  This  is  actually  one  way  of  identifying  pseudotones 
for  invoking  the  provision  in  paragraph  B36.3  of  FAR  Part  36. 

Figure  1-3  also  shows  that  pseudotones  appearing  in  real  data  are  never 
as  pronounced  as  indicated  by  theoretical  simplified  considerations.  Sharp  spectral 
dips  are  "washed  out"  due  to:  (I)  atmospheric  turbulence  (randomizes  propagation 
speeds  and  therefore  the  effective  path  length  difference  or  the  time  delay  of  the 
reflected  relative  to  the  direct  wave),  (2)  surface  impedance  (may  not  be  constant 
as  the  area  of  the  surface  that  reflects  the  sound  towards  the  microphone  moves 
with  the  position  of  the  aircraft),  (3)  surface  randomness  (surface  is  never  ideally 
flat),  and  (4)  the  finite  bandwidth  of  the  analysis  filter. 

Pseudotones  may  have  two  effects  on  the  Effective  Perceived  Noise 
Level  (EPNL)  calculated  from  the  series  of  spectra  of  a  flyover.  One,  the 
Perceived  Noise  Level  (PNL)  at  any  instant  may  be  different  from  that  of  free 
field  conditions.  Two,  the  tone  correction  to  be  added  to  the  PNL  may  be  different 
due  to  the  presence  or  absence  of  pseudotones.  Concerning  the  latter,  FAR  Part 
36  allows  ignoring  the  portion  of  the  spectrum  containing  pseudotones  under 
certain  circumstances  (see  Section  1.2).  Nowadays  it  is  common  practice 
(approved  by  FA  A)  to  ignore  tone  corrections  below  800  Hz,  but  SPL  values  must 
be  included  in  the  noy  calculation.  This  procedure  has  arisen  from  practical 
necessity,  the  800  Hz  cutoff  being  arbitrary.  As  shown  in  Figure  1-3,  pseudotones 
at  higher  frequencies  can  occur  quite  easily. 

Subjecting  the  entire  spectrum  to  tone  correction  brings  up  another 
problem:  the  tone  correction  procedure  prescribed  by  FAR  Part  36  is  considered  by 
many  not  as  advanced  as  the  procedure  recommended  by  SAE  ARP  107 1.^  The 
former  represents  a  stage  in  the  development  of  ARP  1071  and  employs  a  10-step 
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procedure,  whereas  SAE  ARP  1071  employs  a  7-step  procedure  which  gives  more 
accurate  results  in  some  cases  which  are  of  importance  when  considering  pseudo¬ 
tones.  both  methods  usually  give  very  similar  tone  corrections,  except  when  a  tone 
is  shared  by  contiguous  frequency  bands:  both  tone  correction  procedure  require 
that  the  given  spectrum  be  smoothed  so  that  the  difference  due  to  the  tone  can  be 
determined.  The  10-step  procedure  does  not  result  in  a  spectrum  smoothed  as  well 
as  by  the  7-step  method.  This  results  in  a  smaller  tone  correction  for  the  10-step 
method.  Considering  that  pseudotone  frequencies  continually  vary  during  an 
aircraft  flyover  (Figure  1-3),  cases  when  tones  are  shared  by  contiguous  bands  will 
almost  certainly  occur.  Therefore,  if  the  7-step  procedure  were  incorporated  into 
FAR  Part  36  in  place  of  the  10-step  procedure,  the  pseudotone  problem  could  be 
accentuated  because  the  7-step  procedure  recognizes  "tones"  more  readily. 

Another  influence  of  pseudotones  on  EPNL  is  as  follows:  Some  aircraft 
noise  spectra  contain  real  tones  (for  example,  buzz-saw  noise  from  supersonic  tip 
speed  fan  engines  without  inlet  guide  vanes  IGV)  in  the  frequency  region  where 
pseudotones  are  strongest.  These  would  be  ignored  if  the  800  Hz  cutoff  is  blindly 
applied.  They  can  be  distorted  by  pseudotones  and  inaccurate  pure  tone  correc¬ 
tions  applied. 

It  could  be  argued  that  pseudotones  are  really  heard  by  a  person  standing 
or  sitting  outdoors  and  that,  therefore,  they  should  not  be  eliminated  from 
measured  spectra.  However,  in  urban  and  suburban  areas,  people  spend  most  of 
their  time  indoors.**  The  response  of  buildings  to  acoustical  excitation  generally 
favors  the  lower  frequencies,  the  region  where  pseudotones  are  found.  Therefore, 
accuracy  of  the  spectrum  of  the  sound  as  seen  by  the  building  is  desirable 
particularly  in  the  lower  regime.  Inside  a  building,  pseudotones  may  also  occur  but 
will  be  of  a  totally  different  nature  due  to  the  presence  of  so  many  internal 
reflecting  surfaces  in  close  proximity.  In  any  event,  such  effects  are  not  pertinent 
to  aircraft  noise  certification  under  current  rules. 

Without  specifying  in  detail  the  reflecting  surface  for  a  certification 
measurement,  the  applicant  for  an  aircraft  noise  type  certificate  is  free  to  choose 
a  fairly  absorbing  surface  such  as  lush  grass  land  (FAR  Part  36's  stipulation  of  "no 
excessive  sound  absorption  characteristics"  is  not  quantitative  and  therefore  not 
enforceable)  in  order  to  minimize  the  influence  of  the  reflected  wave  resulting  in 
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slightly  lower  noise  levels  than  with  a  harder  surface.  Another  applicant  may  not 
be  able  to  locate  such  favorable  conditions  so  that  aircraft  would  be  compared 
under  different  reference  conditions  which  is  unfair.  Many  of  the  subsequently 
discussed  methods  of  dealing  with  pseudotones  automatically  carry  with  them  the 
benefit  of  removing  this  inequity. 

1.2  Current  FAR  Part  36  Requirements  on  Ground  Reflections 

29 

Applicable  excerpts  from  the  current  version  of  FAR  Part  36  are: 

o  Paragraph  A36. 1  (b): 

"Locations  for  measuring  noise  from  an  aircraft  in  flight  must  be 
surrounded  by  relatively  flat  terrain  having  no  excessive  sound 
absorption  characteristics  such  as  might  be  caused  by  thick, 
matted,  or  tall  grass,  shrubs,  or  wooded  areas." 

o  Paragraph  A36.3(f): 

"The  microphones  must  be  placed  so  that  their  sensing  elements 
are  approximately  4  feet  above  ground." 

o  Paragraph  B36.5  (tone  correction  procedure): 

"For  any  i-th  one-third  octave  band,  at  any  k-th  increment  of  time, 
for  which  the  tone  correction  factor  is  suspected  to  result  from 
something  other  than  (or  in  addition  to)  an  actual  tone  (or  any 
spectral  irregularity  other  than  aircraft  noise),  an  additional  analy¬ 
sis  may  be  made  using  a  filter  with  a  bandwidth  narrower  than  one- 
third  of  an  octave.  If  the  narrow  band  analysis  corroborates  that 
suspicion,  then  a  revised  value  for  the  sound  pressure  level,  SPL"  (i, 
k),  may  be  determined  from  the  analysis  and  used  to  compute  a 
revised  tone  correction  factor,  F(i,  k),  for  that  particular  one-third 
octave  band." 

From  these  inclusions,  it  is  seen  that  FAR  Part  36  recognizes  spectral 
irregularities  other  than  those  from  the  source  (the  aircraft),  allows  that  they  be 
ignored  for  calculating  tone  corrections  to  the  Perceived  Noise  Level,  but 
otherwise  does  not  require  any  procedures  that  would  ovoid  or  compensate  for 
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pseudotones.  It  is  probably  safe  to  soy  that  pseudotones  are  present  to  some  extent 
in  all  certification  tests  because  FAR  Part  36  requires  that  the  microphone  be  1.2 
m  (4  feet)  above  a  ground  that  is  not  overly  absorptive,  i.e.,  quite  capable  of 
reflecting  sound  and  thereby  generating  pseudotones.  The  nature  of  the  ground 
cover  need  only  be  reported  qualitatively. 

In  view  of  the  many  other  corrections  to  the  measured  data  required  by 
FAR  Part  36  (airplane  position,  weather,  ambient  noise)  which  are  generally 
applied  with  a  relatively  high  degree  of  accuracy,  it  appears  very  desirable  to  also 
minimize  the  influence  of  the  ground  surface  and  correct  to  free  field  conditions 
(i.e.,  simulating  the  absence  of  the  ground)  because  the  nature  of  the  ground 
surface  strongly  determines  the  magnitude  of  its  influence.  Particularly  in  the  low 
frequency  regime,  ground  reflections  may  distort  spectra  to  a  much  larger  degree 
than  other  factors  (weather,  ambient  noise,  airplane  position).  An  option,  not 
considered  further  in  this  report,  would  restore  a  +3  dB  correction  to  "free  field" 
data  to  achieve  a  level  which  accounts  only  for  a  uniform  energy  addition  of  the 
ground  reflected  wave  at  all  frequencies. 

1.3  Analytical  Corrections 

We  define  analytical  correction  procedures  as  those  procedures  which 
attempt  to  reconstruct  free  field  flyover  spectra  from  measured  spectra  contami¬ 
nated  by  surface  reflection. 

Table  I- 1  lists  the  analytical  pseudotone  correction  techniques  that  were 
identified  as  candidates  for  thorough  investigation.  Each  technique  is  briefly 
discussed  in  the  following  paragraphs. 

I.  Spectrum  Smoothing 


This  technique  would  apply  a  calculation  procedure  to  the  lower  fre¬ 
quency  portion  of  the  spectrum  that  would  even  out  the  peaks  and 
troughs  caused  by  the  pseudotones.  One  way  would  be  to  assume  that 
the  lower  portion  contains  only  jet  noise  for  which  the  spectrum  shape  is 
very  well  known.  The  spectrum  could  then  be  extrapolated  on  that  basis. 
Another  way  would  be  to  postulate  a  simple  polynomial  expression  for 
the  lower  spectrum  which  could  then  be  least  squares  fitted.  The  danger 
would  of  course  be  that  real  tones  emitted  by  the  aircraft  would  also  be 
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smoothed  over.  Thus,  this  technique  could  probably  only  be  used  reliably 
for  jet  and  fan  jet  aircraft  which  do  not  exhibit  strong  spectral 
variations  for  frequencies  where  pseudotones  occur  (below  I  kHz  or  so), 
such  as  is  the  case  with  multiple  pure  tones  generated  by  transonic  tip 
speed  fans  on  some  older  wide  body  aircraft.  In  particular,  the  technique 
could  not  be  applied  to  general  aviation  propeller  and  rotary  wing 
aircraft  which  emit  most  of  their  acoustic  energy  in  the  frequency 
region  strongly  influenced  by  pseudotones. 

2.  Frequency  Cutoff  for  Tone  Correction 

This  is  the  technique  now  sanctioned  by  FAR  Part  36  in  Appendix  B, 
Paragraph  B36.5(m)  (see  Section  1.2  in  this  report).  Below  a  certain 
frequency  (800  Hz),  any  spectral  peaks  giving  rise  to  a  tone  correction  to 
PNL  may  be  ignored  in  the  calculation  of  the  tone  correction.  This 
provision  arose  out  of  practical  necessities  in  aircraft  noise  certifica¬ 
tion.  The  rigorous  theoretician  might  label  this  technique  a  "band-aid 
approach"  since  it  only  avoids  a  problem  rather  than  solves  or  eliminates 
it.  Also,  the  calculation  of  PNL  is  still  based  on  the  contaminated 
spectra.  Nevertheless,  the  method  has  proven  useful  in  practice  so  that 
it  would  deserve  more  detailed  investigation.  The  cut-off  technique 
actually  departs  from  the  definition  of  an  "analytical"  correction  tech¬ 
nique  given  at  the  beginning  of  this  section  since  a  correction  to  free 
field  is  not  attempted. 

3.  Specular  Reflection  Theory 

The  theory  that  describes  the  superposition  of  the  direct  and  obliquely 
reflected  wave  is  invoked  to  calculate  correction  coefficients  to  trans¬ 
form  the  contaminated  spectra  to  the  spectra  corrected  to  free  field. 
The  value  of  this  technique  depends  on  how  faithfully  the  theory  models 
the  real  world.  Provided  the  modeling  is  faithful  enough,  this  approach 
is  very  defensible,  although  the  data  and  computational  requirements  are 
rather  significant.  To  a  large  extent,  the  success  of  the  technique 
depends  on  the  accurate  knowledge  of  the  acoustic  surface  impedance. 
Although  substantial  progress  has  been  made  in  this  area,^’ 
it  remains  to  be  shown  that,  except  for  certain  special  conditions  such  as 
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a  very  hard  surface,  the  impedance  of  any  practical  surface  can  be 
adequately  specified  for  application  to  "approved"  noise  certification 
procedures. 

4.  Cepstral  Techniques 

The  measured  ("contaminated")  spectrum  in  terms  of  SPL  (in  dti)  versus 

frequency  is  subjected  to  a  Fourier  transformation  from  frequency  space 

into  quefrency  space  (quefrency  =  echo  delay  time).  In  quefrency  space, 

2 1 

the  "spectral"  function  is  called  the  "cepstrum."  If  the  spectrum 
subjected  to  the  Fourier  transform  was  contaminated  by  surface  reflec¬ 
tions,  the  cepstrum  will  show  peaks  at  the  quefrencies  corresponding  to 
the  echo  delay  time  of  the  reflected  signal,  and  to  multiples  of  that 
delay.  In  other  words,  the  influence  of  reflections  may  be  more  readily 
recognized  in  the  cepstrum  than  in  the  spectrum  whence  it  is  more  easily 
removed. 


Figure  1-4  depicts  the  process  described  in  the  previous  paragraph 
applied  to  the  reflection  correction  only.  The  cepstrum  was  obtained  by 
numerical  discrete  Fourier  transformation  of  the  finite  record  shown  at 
the  top  of  the  figure  (128  points  on  the  frequency  axis).  Sharp  peaks  (or 
valleys)  do  indeed  occur  in  the  cepstrum  at  the  integer  multiples  of  the 
true  echo  delay  time.  In  the  method  described  in  Reference  21,  the 
cepstrum  of  the  measured  spectrum  is  subjected  to  "comb"  filtering,  i.e., 
the  cepstral  values  at  the  echo  delay  time  quefrency  (and/or  in  its 
vicinity)  are  set  to  zero.  Subsequent  inverse  transformation  to  the 
frequency  domain  results  in  a  "smoothed"  spectrum. 

Setting  certain  cepstral  values  to  zero  introduces  an  error  in  the 
cepstrum  (and  the  spectrum  after  inverse  transformation)  because, 
strictly  speaking,  only  the  cepstrum  of  the  reflection  should  have  been 
removed,  and  not  also  the  cepstral  values  contributed  by  the  direct  (free 
field)  sound.  If  the  bands  are  spaced  closely  enough,  the  error  is  not 
expected  to  be  significant  for  broad  band  noise  sources.  However,  many 
noise  sources  (propeller,  rotary  wings,  multiple  pure  tones  from  transonic 
fans)  have  tones  which  may  give  rise  to  cepstra  which  look  quite  similar 
to  the  reflection  cepstrum.  During  the  processing  of  0.5-second  flyover 
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data,  the  probability  is  high  that  the  two  cepstra  (reflection  and  source) 
will  or  almost  will  coincide  so  that  comb  filtering  will  remove  the 
greater  part  of  the  acoustic  energy  in  the  signal. 

This  phenomenon  could  be  corrected  by  removing  only  the  reflection 
cepstrum.  However,  the  advantage  of  using  the  simple  comb  filter  is 
lost  and  the  amount  of  computing  work  is  very  similar  to  (or  even 
greater  than)  calculating  the  reflection  correction  directly. 

5.  Combination  Techniques 

These  consist  of  some  combination  of  theoretical  and  empirical  tech¬ 
niques.  Where  theory  is  unable  to  provide  a  realistic  enough  mathe¬ 
matical  model,  empiricism  is  employed  to  complement  the  modeling 
process. 

One  combination  technique  has  been  developed  by  SNECMA  (Appendix  A 
of  Reference  22)  which  examines  each  spectrum  in  the  flyover  series  and 
identifies  the  nulls  (troughs)  presumed  to  be  due  +o  pseudotone  phe¬ 
nomena.  The  location  of  these  on  the  frequency  axis  is  used  to  "tune"  a 
theoretical  reflection  correction  curve.  As  already  pointed  out  for  other 
techniques,  this  technique  is  not  applicable  to  spectra  containing  dis¬ 
crete  frequency  components  in  the  pseudotone  area  (propellers,  rotors, 
transonic  fans)  as  they  have  their  own  peaks  and  troughs.  On  a  one-third 
octave  band  spectrum,  these  true  spectral  peaks  and  valleys  are  difficult 
to  distinguish  even  by  the  trained  eye,  or  require  much  more  sophis¬ 
ticated  analysis  of  the  data  in  the  time  domain  (i.e.,  autocorrelation, 
etc.). 

To  some  extent,  many  so-called  theoretical  methods  incorporate  a 
certain  measure  of  empiricism.  So  does  the  reflection  correction 
technique  described  later  in  Section  1.3.2. 1  (Step  6)  where  the  atmos¬ 
pheric  and  surface  inhomogeneities  and  irregularities  are  accounted  for 
by  a  correlation  coefficient  which  depends  on  an  empirically  determined 
constant. 

1.3.1  Selection  of  Techniques 

For  this  study,  the  following  criteria  were  applied  during  selection  of 
pseudotone  correction  techniques: 
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o  Technical  feasibility  and  accuracy, 

o  Economic  reasonableness, 

o  Enforceability(i.e.,  suitable  for  "approval"  by  regulatory  agencies), 
o  Time  frame  for  implementation. 

The  direct  spectrum  smoothing  technique  is  technically  feasible,  but  its 
accuracy  is  questionable  because  it  would  smooth  over  source  generated  tones.  It 
is  economically  reasonable;  its  application  would  be  straightforward  and  therefore 
cheap.  It  should  be  enforceable  if  a  standard  technique  were  adopted.  It  is 
immediately  implementable. 

The  frequency  cutoff  for  the  tone  correction  technique  has  a  demon¬ 
strated  technical  feasibility.  The  criterion  of  accuracy  is  not  directly  applicable 
because  the  reflections  are  not  removed  (errors  are  indirectly  evaluated  in  Section 
1.3.4).  The  technique's  application  is  very  easy  and  therefore  economically 
reasonable.  It  is  enforceable  if  a  certain  cutoff  frequency  (800  Hz  in  the  current 
FAR  Part  36)  is  chosen;  if  it  varies  (say,  depending  on  elevation  angle),  enforce¬ 
ability  may  suffer.  In  any  event,  this  technique  is  already  implemented  by  FAA. 

The  theoretical  reflection  correction  technique  is  technically  feasible. 
Its  accuracy  is  great  if  the  assumptions  made  in  the  derivation  of  the  theoretical 
formulas  are  reasonably  well  approached.  Concerning  the  assumptions  of  a 
perfectly  still  and  homogeneous  atmosphere  and  of  a  completely  plane  surface,  this 
assumption  is  practically  never  satisfied.  This  deficiency  can  be  compensated  for 
by  an  empirical  method.  The  required  calculations  for  this  theoretical  technique 
are  not  straightforward  but  would  not  add  an  undue  burden  to  certification  noise 
processing.  One-time  development  and  testing  of  the  necessary  computer  program 
code  should  take,  at  the  most,  several  man-weeks. 

The  cepstral  technique  may  be  technically  feasible  but  difficult,  practi¬ 
cally,  since  it  requires  narrow  band  analysis.  Whereas  the  currently  implemented 
one-third  octave  band  analysis  requires  obtaining  24  values  every  0.5  second, 
narrow  band  analysis  would  require  at  least  five  times  as  much  if  not  more.  If 
narrow  band  (i.e.,  fixed  bandwidth  filter)  analysis  was  not  used,  the  data  points 
would  not  be  equispaced  in  the  frequency  domain  and  Fast  Fourier  Transform 
techniques  could  not  be  applied  which  would  make  the  process  more  time- 
consuming  and  therefore  more  expensive.  Accuracy  is  acceptable  for  smooth 
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source  spectra  (such  as  jet  noise),  but  may  not  be  if  tones  with  harmonics  are 
present.  Because  of  the  problems  with  technical  feasibility,  the  technique  would 
probably  be  expensive  in  its  implementation.  An  acceptable,  codified  and 
enforceable  cepstral  pseudotone  correction  technique  would  probably  take  several 
years  to  implement. 

SNECMA's  combined  theoretical  and  empirical  technique  is  not  always 
technically  feasible  as  it  breaks  down  in  the  presence  of  spectral  irregularities  due 
to  the  source.  If  it  is  applicable,  it  is  probably  very  accurate.  Its  economics, 
enforceability  and  time  frame  of  implementation  are  comparable  to  the  theoretical 
reflection  correction  technique. 

As  has  already  been  mentioned  above,  the  latter  may  be  enhanced  by 
empirical  adjustments.  This  technique  shows  the  most  promise  of  providing  a 
rigorous,  accurate,  and  not  too  complicated  pseudotone  correction.  It  is  therefore 
selected  as  the  primary  technique  to  be  investigated  in  this  study. 

As  a  secondary  technique,  the  frequency  cutoff  for  tone  correction  is 
selected  in  order  to  investigate  the  relative  accuracy  of  an  already  codified 
provision. 

The  other  techniques  (direct  spectrum  smoothing,  cepstrum,  SNECMA) 
are  rejected  for  the  purposes  of  this  study  primarily  because  of  their  inability  to 
properly  deal  with  spectral  irregularities  emanating  from  the  source  which 
superimpose  themselves  over  the  pseudotones. 

1.3.2  Ground  Reflection  Correction 

1.3.2. 1  Theory 

Recently,  an  easy-to-use  and  effective,  yet  theoretically  well  founded 
procedure  has  been  published  for  calculating  ground  reflection  effects.'*’*’  This 
method  is  adopted  here  and  briefly  summarized  in  the  following  paragraphs.  The 
theory  is  based  on  the  following  assumptions: 

o  Ground  surface  is  flat  with  a  clearly  identifiable  interface  between 
air  and  ground  (over  a  lush  grassy  surface  it  is  often  unclear  where 
this  boundary  is). 

o  Ground  surface  is  "locally  reacting"  (i.e.,  sound  waves  transmitted 
into  the  ground  propagate  only  perpendicularly  to  the  surface). 

o  Atmosphere  is  uniform  (no  wind,  no  temperature  gradients). 


o  Spectrum  levels  vary  slowly  over  any  one  third-octave  band  (this 
condition  may  well  be  violated  in  practice  where  single  tones  are 
present  but  very  narrow  band  analysis  would  be  required  to  resolve 
this  problem). 

o  Ground  impedance  is  large  relative  to  that  for  air  (more  details  on 
this  assumption  are  given  later). 

The  geometry  of  the  situation  to  be  modelled  is  shown  in  Figure  I  -  I. 
The  calculation  framework  for  arriving  at  the  ground  reflection  correction  is  given 
in  the  following  steps.  This  correction  must  be  added  to  a  measured  signal  in  order 
to  remove  the  ground  reflection  effect. 

Calculation  Framework: 

Step  I:  Determine  Ground  Complex  Impedance  or  Admittance  =  I /impedance  as 

a  function  of  frequency:  A  separate  section  (1. 3.2.2)  is  devoted  to  this 
subject. 

Step  2:  Determine  Ambient  Speed  of  Sound  Cq. 

a.  Approximate  Method: 

Cq  =  331.6  (I  +  tc/273.l5)^  in  meters/second, 
where  tc  is  the  temperature  in  degrees  Celsius. 

b.  Thermodynamic  Method  (includes  effect  of  humidity): 

The  method  employed  here  for  accounting  for  effects  of  humidity 
on  the  ambient  speed  of  sound  is  chosen  for  computational  conveni¬ 
ence.  An  alternative,  semiempirical  method,  based  on  measured 
data,  is  also  available  in  Appendix  A  of  Reference  30.  The  two 
methods  differ  by  less  than  0.1%. 

t^  =  t  +  273.15,  temperature  in  degrees  Kelvin; 

tf  =  t^/273.16,  triple  point  ratio; 

Determine  water  vapor  partial  pressure  divided  by  total  pressure: 

\/p  -  __!  m  exponent 

vr  "  100  p  IU 
o 

where:  RH  =  relative  humidity  in  percent, 
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kji  if 


characters  denote  vector  quantities),  the  unit  vector  u  in  direction  of 
the  flightpath,  and  aircraft  speed  S.  With  time  denoted  by  t,  tQ 
indicating  the  time  at  the  apparent  location,  points  on  the  flightpath 
are  given  by  A  +  uS  (t  -  t  ).  With  t  denoting  the  time  at  emission  the 
following  relationship  holds  for  the  emission  location  E  (see  Figure 
1-4.5): 

E=  A+uS(te-tQ) 

Also,  t  -  tg  must  equal  the  sound  travel  time  along  E : 


which  may  be  substituted  into  the  previous  equation  yielding: 

A  =  E  +  o  M  J  E  J 

where  M  is  the  flight  Mach  Number  S/co. 

This  equation  may  be  solved  for  E.  The  solution  is  obtained  as  follows 
(for  the  derivation  see  Appendix  B): 

Given:  A  =  (A j ,  A2,  A3,),  M,  u  =  (uj,  U2,  u3,)  (I-  and  2-direction  in 
ground  plane,  3-direction  vertical) 

Calculate: 


B2  —  A2  —  A  |  j ,  =  A3  *"  ^  |  | 

2  2 
C  |  =  A  |  +  u  |  U2  82  M  +  u  |  u3  B3M 

C2  =  A^-Uj2  M2  (B22  +  B32) 

E,  =(C,  -sign(u()(c,2-C2  (I  -  M2))fe)/(l  -  M2) 
E2  =  B2  +  E|  u2/uj 

E3  =  B3  +  EI  u3/uI 


where 


sign  (u  | )  = 


I  for  u  j  >0 
-I  for  u |  <  0 


The  case  of  Uj  =  0  (or  almost  zero)  needs  special  consideration  because 
it  is  a  devisor  in  the  above  calculations,  and  an  argument  of  the  sign 


Figure  1-4.5.  Illustration  of  Apparent  and  Sound  Emission  Aircraft  Location 
Vectors 


pQ  =  atmospheric  pressure  in  units  of  atmospheres 
(I  atmosphere  =  1.01325  x  10^  N/m^) 
exponent  =  8.58  -  (l0.8/tr)  -  5.03  log  ,Q  (tr)  +  23  x  IO“4,77/tr 
Determine  moisture  content  X  in  kg  of  water  per  kg  of  dry  air: 

X  =  0.6220  VP/(I-VP) 

where  0.6220  is  the  ratio  of  the  air  and  water  vapor  gas  constants. 
Determine  the  gas  constant  of  the  oir-water  vapor  mixture: 

Rm  =  (47.06  X  +  29.27)/(l  +  X) 

where  47.06  is  the  water  vapor  gas  constant,  and  29.27  is  the  air 
gas  constant,  both  in  units  of  meters/degrees  Celsius. 

Determine  the  ratio  of  specific  heats  k  for  the  vapor-  air  mixture: 

k  =  1.402  (X  +  0.868)/(  1 .0526  X  +  0.868) 

where  0.868  is  the  air-to-water  ratio  of  specific  heat  at  constant 
pressure,  and  1.0526  is  the  air-to-water  ratio  of  the  ratios  of 
specific  heats  (=  1.402/1.332). 

Determine  the  ambient  speed  of  sound  cQ: 

LL 

c0  =  (k  9  in  meters/second, 

2 

where  g  =  9.81  m/sec  =  gravitational  acceleration. 

Determine  the  Aircraft  Emission  Location.  The  time  when  a  sound 
emitted  by  the  aircraft  is  heard  (apparent  location)  and  the  time  when 
this  sound  was  emitted  differ  by  the  sound  travel  time  from  the 
emission  location  to  the  observer.  Almost  all  of  the  time,  we  are  given 
the  apparent  location  and  the  aircraft  speed  so  that  the  emission 
location  must  be  determined  in  order  to  know  the  correct  angle  of 
sound  incidence. 


Without  loss  of  generality,  the  origin  of  the  coordinate  system  may  be 
assumed  to  coincide  with  the  observer.  The  flightpath  is  assumed  to  be 
a  straight  line  specified  by  the  apparent  aircraft  location  A  (bold 


function  which  is  undefined  if  the  argument  equals  zero.  Since  we  can 
safely  exclude  the  case  U|  =  u2  =  6  (vertical  flight  path)  it  is  only 
necessary  to  rotate  the  coordinate  system  around  the  3-direction  by  90 
degrees  for  the  duration  of  the  emission  location  calculation. 

Knowing  the  components  of  E,  sound  travel  distance  R  and  elevation 
angle  b  are  easily  calculated: 

R  =  (E12  +  E22  +  E32)fe 

b  =  sin"1  (E3/R) 

The  angle  of  incidence  0  is: 

0  =  90°  -  b 

Determine  the  Path  Length  Difference. 

The  path  length  difference  d  between  direct  and  reflected  wave  is  r2  - 
r  | .  However,  since  r2  =  r  | ,  which  is  much  greater  than  h,  d  is  more 
conveniently  calculated  from 

d  =  2  h  cos  0 


Determine  Complex  Reflection  Coefficient.  Define  a  complex  variable 
t  as  follows: 

t  =  (cos  0  +  n)  (k  R/2i)k  ,  i  =  +  (-1)^ 

where:  n  =  complex  admittance  (from  Step  I) 

k  =  wavenumber  =  2  TT  f/cQ  (  TT  =  3.141592...  ) 
f  =  frequency 

Then  calculate 

F(t)  =  I  -  t  W  (it) 


where  W  is  the  complex  error  function.  Appendix  C  gives  a  detailed 
calculation  procedure  for  the  function  F(t),  including  a  computer 
program. 


The  plane  wave  complex  reflection  coefficient  is: 

Gp  =  (cos  0  -  n)/(cos  0  +  n) 


f 


The  general  complex  reflection  coefficient  is: 

Gc' V('-Gpm”'-G'ia 

where:  G  =  I  G  I  ,  a  =  tan“ 1  (Im  G  /Re  G  ) 

I  C  |  C  C 

Determine  Correlation  Coefficient  Adjustment  to  account  for  Band 
Averaging  of  Finite  Bandwidths  of  Noise 

The  band  averaging  factor  B  is  given  by: 

B  =  (sin  (c  k  d/2))/(fckd/2) 

where  c  is  the  filter  bandwidth  relative  to  the  band  (geometric)  center 
frequency.  For  one-third  octave  bands  (using  the  approximate  equiva¬ 
lent  of  one-tenth  decade  bands): 


c  =  10  0,05  -10  ‘0,°5  =  0.23077 


Determine  Correlation  Coefficient  (correction  for  non-ideal  atmos¬ 
phere). 


Inhomogeneities  in  the  atmosphere  and  the  reflecting  surface  cause 
amplitude  and  phase  differences  between  direct  and  reflected  signals  to 
vary  randomly.  A  correction  factor  to  the  correlation  coefficient  is 
used  to  account  for  this  phenomenon  empirically.  Higher  frequencies 
(wave  numbers)  will  be  affected  more  strongly  than  lower  ones  and 
coherence  between  the  signals  should  also  be  expected  to  decrease  with 
path  length  difference  d.  Reference  6  suggests  the  following  form  for 
this  correlation  correction  factor  C  which  is  adopted  here: 


C  =  exp  -(,inc  kd> 


where  f-  is  an  incoherence  factor  with  a  recommended  value  of  0.01 
(subject  to  empirical  adjustment).  Application  of  this  correlation 
coefficient  relaxes  some  of  the  assumptions  made  in  the  beginning  of 
this  section. 


Determine  Reflection  Correction 

The  reflection  correction,  in  decibels,  that  must  be  added  to  a  band 
level  in  order  to  remove  the  ground  reflection  effect  becomes: 
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pc  =  -  »0  log10  (l  ♦  G2  ♦  2  G  B  C  cos  (o  ♦  k  d)) 

An  additional  correction  factor  should,  theoretically,  be  included  in  the 
argument  of  the  cosine  term  in  this  expression;  however,  it  is  essentially  equal  to 
unity  for  one-third  octave  bands  and  is  neglected  here.  The  effects  of  source 
correction  due  to  the  finite  averaging  of  the  source  spectra  have  been  neglected 
here.  It  has  been  shown  that  averaging  times  of  0.5  seconds  produce  a  difference 
of  only  0.2  dB  by  including  this  extra  term,  *7at  least  at  high  angles  of  incidence. 

1. 3.2.2  Complex  Ground  Impedance 

The  accurate  prediction  of  outdoor  ground  impedance  on  the  basis  of 
measured  data  and/or  theoretical  models  is  generally  difficult  due  to  the  hetero¬ 
geneous,  anisotropic  nature  of  most  surfaces.  Impedance  (or  its  reciprocal, 
admittance)  is  a  complex  quantity  generally  varying  as  a  function  of  both 
frequency  and  angle  of  incidence.  Numerous  studies  have  attempted  to  evaluate 
the  effect  of  finite  ground  impedance  indirectly  by  experimental  or  analytical 
studies  of  horizontal  sound  propagation  or  have  measured  impedances  for  various 

surfaces  directly  or  using  an  assortment  of  elaborate  measurement  tech- 
I  I  19 

niques.  Invariably,  certain  assumptions  must  be  made  in  order  to  simplify  the 
modeling  process:  all  surfaces  are  assumed  to  be  homogeneous,  isotropic,  and 
"locally  reacting." 

12 

Delany  and  Bazley  were  the  first  to  define  normalized  impedance 
(normalized  with  respect  to  air)  as  a  function  of  frequency,  f,  and  flow  resistivity, 
s,  for  fibrous  absorbent  materials  (cotton,  fiberglass,  etc.).  Both  of  these 
quantities  are  easily  obtained  for  many  materials  of  interest.  Indeed,  the  power 
law  relations  developed  by  Delany  and  Bazley: 

R/o  cQ  =  I  +  9.08  (f/s)-75 

X//>oco=  1 1  -9  (f/s)"*73  * 

agree  well  with  the  measured  data  for  the  materials  tested.  Here  R  represents 
the  real  part  and  X  the  imaginary  part  of  the  surface  impedance,  and  yoQ  cQ  the 
specific  impedance  of  air  at  standard  atmospheric  conditions. 


In  the  present  formulation  of  the  theory,  the  imaginary  part  of  the  impedance  must 
have  a  positive  sign.  Some  references  give  it  as  negative.  See  Reference  31  for  a 
discussion. 
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More  recently,  other  studies  have  shown  that  these  power  law  relations 
agree  fairly  well  with  measurements  taken  over  grass, hard  ground,*^  spaded 
soil,  and  snow.  Although  measured  and  extrapolated  values  of  impedance  were 
found  to  be  in  excellent  agreement  with  the  above  model  for  values  of  (f/s) 
between  0.5  and  1000  cgs  rayls/cmj'*  caution  should  be  used  when  applying  the 
model  to  any  surface.  Moisture  content,  root  structure,  and  the  geological 
characteristics  of  the  subsoil  (and  topsoil)  play  an  important  role  jn  describing  the 
flow  resistivity  of  the  surface  material.  For  values  of  s  for  common  surface 
materials,  see  Refererence  20. 
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1.3.2. 3  Application  of  Ground  Reflection 

A  computer  program  was  written  to  apply  the  calculation  of  ground 
impedance  and  reflection  correction  described  in  the  previous  sections.  In  order  to 
demonstrate  the  dependance  of  the  reflection  correction  on  ground  impedance, 
Figures  1-5  through  1-8  were  created  by  the  computer  program.  In  all  four  cases, 
the  reflection  correction  (i.e.,  the  amount  in  dB  that  is  to  be  added  to  measured 
spectra  in  order  to  remove  the  ground  effect)  is  plotted  as  a  function  of  frequency 
(a  point  at  each  1/3  -  octave  band  center  frequency)  and  time  (a  spectrum  every  Ifc 
seconds).  An  aircraft  is  assumed  to  execute  a  level  flyover  (altitude  1000  feet 
(304.8  m))  at  a  constant  speed  of  160  knots  (82.3  m/s).  Visual  overhead  occurs  at 
10  seconds  into  the  20  sec  display  window  (actual  or  emission  overhead  occurs 
about  I  second  later).  The  incoherence  factor  (f.  )  was  set  to  0.01,  temperature 

to  20°  Celsius,  humidity  to  70%,  and  pressure  to  I  atmosphere.  The  following  four 
extreme  values  of  normalized  impedance  were  evaluated  to  illustrate  the  general 
influence  of  this  parameter. 


Figure 

Complex  Impedance 

1-5 

1-6 

l-i  | 

1  +  i  \  constant  for  all  frequencies 

1-7 

1000  +  iO  j 

1-8 

frequency  dependent  with  flow  resistivity 
s  =  100  cgs  rayls/cm 

The  normalized  impedances  used  in  Figure  I -8  are  listed  in  Table  1-2. 

All  four  spectral  time  histories  of  the  reflection  correction  exhibit  the 
expected  shift  of  the  peaks-valleys  pattern  towards  lower  frequencies  in  the 
vicinity  of  the  overhead  position.  Figure  1-7  is  most  easily  interpreted  since  it 
represents  very  nearly  the  case  of  a  perfect  reflector.  As  should  be  expected,  the 
reflection  correction  tends  towards  -6  db  for  low  frequencies  (acoustic  pressure 
doubling  near  perfect  reflector),  and  towards  -3  dB  at  high  frequencies  (acoustic 
intensity  doubling).  The  impedance  used  in  Figure  1-5  (l-i)  is  that  of  a  very 
absorbent  and  compliant  surface.  Because  of  the  high  absorption,  pressure  and 
intensity  do  not  increase  nearly  as  much  as  near  a  hard  surface.  Changing  the  sign 
of  the  reactance  (impedance  =  I  +  i,  Figure  1-6  causes  a  radical  change  in  the 
reflection  correction  pattern;  that  is,  however,  not  a  realistic  case  as  practical 
surfaces  have  a  positive  real  part  and  a  negative  (i.e.,  stiffness)  imaginary  part  of 
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CONSTANT  80  08  FOR  REFLECTION  CORRECTION  OISPLBT 


0  dB  for 
Top  Curvi 


Figure  1-7.  Plot  of  Spectral  Time  History  of  Reflection  Correction:  Level  (304.8  m) 
Constant  Speed  (82.3  m/sec)  Flyover,  Visual  Overhead  at  10  Seconds, 
Constant  Normalized  Impedance  =  1000  +  iO. 


CONSTANT  80  08  FOR  REFLECTION  CORRECT 
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FREQUENCY.  HZ 
2  3  «  5  6  7  B 
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3  3  «  5  8  7  8  9  A  O’ 
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Figure  1-8.  Plot  of  Spectral  Time  History  of  Reflection  Correction:  Level  (304.8  m) 
Constant  Speed  (82.3  m/sec)  Flyover,  Visual  Overhead  at  10  Seconds, 
Normalized  Impedance  Depends  on  Frequency  (Calculated  with  Flow 
Resistivity  =  100  cgs  rayls/cm).  Imaginary  Part  of  Ground  Impedance 
is  Negative.  See  Footnote  on  Page  1-22. 
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Normalized  Impedance  and  Admittance  as  a  Function  of  Frequency  as  Used  in 
Figure  1-8,  with  Ground  Surface  Flow  Resistivity  =  100  cgs  rayls/cm 
Imaginary  Part  of  Complex  Impedance  is  Negative  (See  Footnote  on  Page  1-22) 

(Computed  from  Equation  at  Beginning  of  Section  1 .3.2.2) 
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the  impedance.  By  contrast,  Figure  1-8  presents  the  reflection  correction  for  a 
practical  surface  with  a  flow  resistivity  of  100  cgs  rayls/cm  which  is  probably 
representative  of  lush  grass  over  a  very  porous  and  dry  soil.  These  corrections  tend 
toward  -6  dB  at  low  frequencies  as  should  be  expected  due  to  the  relatively  high 
impedance  there  (Table  1-2).  At  high  frequencies,  the  corrections  tend  toward  0 
(1J  due  to  the  associated  low  impedances.  Comparison  with  Figure  1-7  shows  that 
there  is,  apart  from  the  high  frequency  3  dB  shift,  little  difference  between  the 
corrections  for  a  soft  grassy  surface  and  a  hard  surface  in  the  vicinity  of  the 
overhead  positions  (small  angle  of  incidence).  Significant  differences  appear 
upwards  on  the  figure  (earlier  in  time)  at  the  fifth  curve  from  the  top  (corres¬ 
ponding  to  an  angle  of  incidence  of  55  to  60°)  and  downward  to  the  second  or  third 
curve  from  the  bottom  (corresponding  to  the  same  approximate  angle  of  incidence). 
This  would  indicate  that  it  is  not  vital  to  know  the  exact  surface  impedance,  as 
long  as  some  reasonable  value  is  used,  for  near  overhead  flyovers  since  the  greatest 
contribution  to  single  event  metrics  (such  as  EPNL)  comes  from  the  noise  emitted 
near  overhead.  For  sideline  measurements,  the  angle  of  incidence  (O  in  Figure  l-l) 
will  hardly  ever  be  less  than  30°  (for  a  light  B727)  and  often  not  go  below  55°  (for  a 
heavy  B747).  Accurate  knowledge  of  ground  impedance  then  becomes  important. 

When  we  attempted  to  apply  the  reflection  correction  to  measured  data, 
serious  problems  were  encountered.  Before  these  problems  are  discussed  in  detail, 
the  results  of  several  correction  attempts  are  presented.  Appendix  A  contains 
graphical  representations  of  almost  all  the  measured  aircraft  noise  events  utilized 
in  this  study:  the  spectral  time  histories  between  the  10-dB-down  points  are 
shown.  The  following  flights  were  selected  for  attempting  the  reflection  correc¬ 
tion: 

A  B727  on  approach  (Flight  I  in  Appendix  A) 

A  DC- 10  on  approach  (spectral  time  history  plot  not  in  Appendix  A  but  in 
this  section); 

A  B727  on  takeoff  (Flight  G  in  Appendix  A). 

(a)  B727  on  Approach 

Figure  A-13  (Appendix  A,  page  A-30)  shows  the  spectral  time 
history  for  the  10-dB-down  interval  as  measured  by  a  microphone  1.2  m 
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over  soft  ground.  The  pattern  created  by  reflection  interference  is 
easily  discernible.  The  frequency  shift  of  the  pattern  with  time  can  also 
be  gleaned  from  the  figure.  The  general  position  of  the  aircraft  was  not 
measured  except  that  the  time  of  exact  overhead  was  noted.  It  was  then 
assumed  that  the  aircraft  was  on  the  3°  ILS  glide  slope.  The  altitude  at 
overhead  could  then  be  calculated  knowing  the  distance  of  the  micro¬ 
phone  from  the  runway  threshold.  The  speed  of  the  aircraft  was  not 
measured;  a  nominal  speed  of  72  m/sec  (140  knots)  was  used.  With  an 
incoherence  factor  of  0.1  and  a  surface  flow  resistivity  of  100  cgs 
rayls/cm  the  correction  was  then  carried  out.  The  result  of  Figure  1-9 
shows  that  the  attempt  was  to  a  large  extent  successful.  The  character¬ 
istic  pseudotone  pattern  is  removed.  One  should  expect  that  the 
corrected  STH  will  be  very  similar  to  the  ones  measured  for  the  same 
flight  at  10  m  and  on  the  ground.  Comparison  with  Figures  A-ll  (page  A- 
28)  and  A-15  (page  A-32)  bears  out  this  similarity  except  for  some 
"waves"  in  Figure  1-9  in  the  5th  and  6th  spectrum  from  the  top  in  the 
lower  frequency  portion  which  are  not  present  in  Figure  A-ll  and  A-15. 

(b)  DC- 1 0  on  Approach 

In  this  case,  the  aircraft  flew  on  an  approach  path  for  which  the 
ground  projection's  perpendicular  distance  to  the  measuring  station 
was  230  m  (750  ft).  The  elevation  angle  was  about  30°  (it  is  90° 
for  overhead  flights).  The  measured  STH  with  pseudotones  is 
shown  in  Figure  I -10.  The  attempt  to  remove  pseudotones  is 
shown  in  Figure  I -I  I.  There  is  only  one  spectrum  (the  4th  from 
the  bottom)  for  which  pseudotones  are  almost  removed.  For  the 
rest  of  the  STH,  the  attempt  is  a  failure  for  which  there  are  these 
possible  explanations: 

o  The  exact  time  of  closest  approach  was  not  known  with 
sufficient  accuracy.  This  may  cause  a  shift  in  the  correction 
pattern  and  result  in  the  "out-of-phase"  appearance  of  Figure 
l-ll. 

o  The  aircraft's  speed  was  substantially  different  from  the  one 
assumed  for  the  correction  attempt  (speed  was  not  measured 
during  the  acoustic  data  acquisition). 
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Figure  1-9.  Spectral  Time  History  of  a  B727  Overhead  Approach 
after  Reflection  Correction;  Flow  Resistivity 
100  cgs  rayls/cm,  Incoherence  Factor  0.1 
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Figure  1-10.  Spec  fra  I  Time  History  for  the  10-dB-down  Interval  of  o  DC-10  on  o 
Sideline  Approach,  Microphone  1.2  m  over  Soft  Ground,  Measured 
Data  Smoothed  in  the  Time  Domain.  Dataset  Sentinel  9000010010. 


(c)  B727  on  Takeoff 


Figure  A-G2  (page  A-80)  shows  the  spectral  time  history  for  the 
10-dB-down  interval  as  measured  by  a  microphone  1.2  m  over  hard 
ground.  The  reflection  interference  pattern  is  easily  discernible. 
Aircraft  position  was  obtained  through  the  FAA's  Aircraft  Radar 
Tracking  System  (ARTS)  at  LAX  (Los  Angeles  International  Air¬ 
port).  The  result  of  the  correction  attempt  is  shown  in  Figure 
1-12.  It  may  be  seen  that  there  was  moderate  success  in  the  lower 
(i.e.,  later)  part  of  the  spectral  time  history  in  removing  the 
interference  pattern.  However,  in  the  early  portion  (top  part  of 
Figure  1-12)  the  correction  is,  so  to  speak,  "out  of  phase,"  i.e., 
increasing  the  peaks  and  troughs  rather  than  removing  them.  The 
lower  part's  moderate  success  indicates  that  the  correction's 
magnitude  was  about  right,  and  that  aircraft  position  was  known 
with  sufficient  accuracy  for  the  part  of  the  time  history  when  the 
aircraft  was  receding  from  the  microphone.  The  failure  in  the 
early  part  was  probably  caused  by  inaccuracies  in  aircraft  position 
data:  the  ARTS  radar  provides  position  data  at  intervals  of 

somewhat  less  than  5  seconds;  in  the  critical  area  (top  of  Figure 
1-12)  several  ARTS  points  were  missing  so  that  the  aircraft's  flight 
path  was  not  known  with  sufficient  accuracy.  However,  other 
problems  as  discussed  for  the  DC- 10  approach  may  also  play  a  role. 

Discussion  of  Problems 

In  the  present  form,  the  analytical  reflection  correction  procedure  does 
not  perform  satisfactorily.  While  we  believe  that  there  is  no  fundamental  error  in 
the  analytical  procedures,  they  do  not  take  account  of  several  physical  phenomena 
which  may,  singly  or  together,  cause  significant  departures  from  the  assumed 
idealized  conditions  of  the  analysis: 

o  Refraction  of  sound  waves  due  to  vertical  temperature  and  wind 
gradients  in  the  atmosphere  may  change  (both  increase  or  de¬ 
crease)  the  actual  angle  of  incidence,  which  results  in  a  pathlength 
difference  different  from  that  of  the  idealized  condition;  hence  a 


spectral  plot  reR  consecutive  ttie  intervals 


shift  on  the  frequency  axis  of  the  interference  pattern  might  be 
observed.  This  phenomenon  is  only  effective  at  large  incidence 
angles,  i.e.,  near  grazing  incidences,  so  that  it  would  not  account 
for  any  frequency  shifts  near  the  overhead  position. 

o  The  finite  extent  of  the  sound  source  would  tend  to  even  out  the 
interference  pattern's  peaks  and  troughs  since  sound  arrives  simul¬ 
taneously  from  a  range  of  incidence  angles.  This  might  account  for 
the  often  observed  overcorrection  of  the  pseudotone  pattern.  Note 
that  for  the  case  of  the  DC- 10  on  approach  the  distance  between 
engines  is  about  1/4  of  the  aircraft's  altitude. 

o  Microphone  height:  for  hard  ground  (asphalt,  concrete)  there  is  no 
doubt  about  how  to  measure  microphone  height.  For  soft  ground, 
there  usually  is  no  clearly  defined  boundary.  What  is  the  "effec¬ 
tive"  height  over  grass?  At  its  tips,  or  at  the  roots?  Notice  tha* 
the  B727  takeoff  correction  attempt  (hard  ground)  was  marginally 
more  successful  than  that  for  the  DC- 10  (grass  surface).  However, 
this  is  in  contrast  to  the  success  of  the  B727  correction  (soft 
ground). 

o  Incorporation  of  the  so-called  "surface  wave"  in  the  theoretical 
model  for  ground  reflection  effects.  This  refinement  is  expected 
to  be  significant  at  low  grazing  angles  (i.e.,  sideline  positions). 

By  contrast,  we  believe  that  ground  admittance  need  not  be  known  with 
great  accuracy  because  both  "soft"  and  "hard"  grounds  appear  "hard"  to  the 
relatively  low  frequencies  at  which  pseudotones  commonly  occur.  At  high 
incidence  angles,  however  (i.e.,  low  elevation  angles),  the  interference  pattern 
moves  to  higher  frequencies  whence  more  accurate  knowledge  of  ground  admit¬ 
tance  would  be  required.  This  may  be  another  explanation  for  the  failure  of  the 
DC- 10  "sideline"  approach  correction.  Finally,  it  should  be  pointed  out  that  the 
only  previously  known  attempt  to  analytically  predict  ground  reflection  effects  for 
aircraft  flyover  signals  also  found  that  predictions  were  frequently  not  in  agree¬ 
ment  with  observations  (Reference  17).  The  one  success  (the  B727  approach), 
however,  encourages  further  future  work  on  the  subject.  Refinements  which  will 
lead  to  more  successes  include: 

o  More  accurate  knowledge  of  source  position  with  time: 

o  Finite  source  extent  effects  taken  into  account. 

o  Incorporate  the  effect  of  surface  waves  in  the  theoretical  model. 
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Position  Accuracy  Requirements  for  Ground  Reflection  Correction 


The  analysis  of  the  previous  sections  demonstrates  that  ground  reflection 
strongly  depends  on  the  angle  of  incidence  which  continuously  varies  during  a 
flyover.  To  a  lesser  extent,  the  correction  also  depends  upon  actual  distance  and 
surface  impedance.  In  this  section,  we  examine  the  accuracy  with  which  the  angle 
of  incidence  should  be  known  in  order  to  perform  a  reliable  correction  for  PNLT 
calculation  purposes.  In  practical  terms,  we  examine  the  accuracy  with  which  the 
noise  emission  location  of  the  aircraft  should  be  determined,  i.e.,  without  explicit 
mention  of  the  angle  of  incidence. 

The  formula  for  the  reflection  correction  is  given  in  Step  8  of  Section 

I. 3.2.1: 

Pc  =  -10  log | q  (I  +  G^  +  2GBC  cos  (a  +  kd)) 

Let  p  l~3  any  position  coordinate  of  the  aircraft.  The  error  in  Pc  due  to  an 
error  p  in  determining  p  may  be  expressed  as: 


In  the  worst  case,  the  largest  error  happens  to  occur  at  a  frequency  with 
a  sound  level  that  dominates  the  PNL  calculation.  In  that  case,  the  final  error  e^ 
would  equal  ep  .  However,  we  must  also  consider  the  tone  correction  to  PNL 
which,  again  in  fhe  worst  case,  can  be  as  much  as  one-third  of  the  amount  that  the 
band  level  exceeds  some  average  of  the  neighboring  bands.  An  upper  bound  for  the 
magnitude  of  the  final  PNLT  error  therefore  is  assumed  to  be: 

e  e 
ef  "  3  P 

c 


This  is  used  in  the  subsequent  discussions. 

Introducing  the  incidence  angle  0  as  an  intermediate  variable,  the  error  e 

becomes: 


e 


_d_y 

d  P 


A  p 
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dPydQ  -was  derived  formally  from  the  equations  in  Section  1.3.2. 1  with  one 
simplifying  assumption:  the  complex  function  F(t)  (Step  5)  equals  I.  The  details  of 
these  derivations  are  not  presented  here.  The  resulting  formulas  and  equations 
were  cast  into  a  computer  program  which  determined  dP^dO  as  a  function  of  angle 
0  and  frequency.  For  many  values  of  0  (I  degree  increments  from  1°  to  85°)  the 
maximum  absolute  value  Q  of  dF^/d0  was  determined.  For  a  typical  acoustically 
soft  surface  (surface  flow  resistivity  of  300  cgs  rayls/cm),  this  is  plotted  in  Figure 
1-13,  and  was  fitted  by  least  squares  by  the  following  function: 

Q  _  |0-l .37602  +  0.682675  x  I0_l  0  -  0.1  17608  x  I0‘2  02  +  0.78/710  x  I0"5  03 
(where  0  is  in  degrees). 

The  same  procedure  was  applied  to  an  acoustically  hard  surface  (surface 
flow  sensitivity  of  ID3  cgs  rayls/cm)  resulting  in  the  following  fit: 

Q.  ,0- 1.4 1 795  +  0.731505  x  I0"1  9  -0.1 35248  x  1 0"2  02  +  0.98001 1  x  I0"5  03 

Rearranging  the  above  equation  for  e  and  choosing  particular  values  for 
a  tolerable  e,  we  may  write 


6  P  = 


Q 


3 

4 


e 


Choosing  the  position  variable  p  to  indicate  either  vertical  or  horizontal  directions, 
graphs  like  the  ones  shown  in  Figures  1-14  and  1-15  result.  The  data  in  these 
figures  were  actually  obtained  directly  from  the  computer  listings  from  which  the 
above  Q-0  relations  were  derived.  The  arrangement  on  these  graphs  was  chosen 
such  as  to  be  able  to  choose  any  maximum  error  e  (numbers  in  dB  on  top  of  middle 
section)  and  immediately  find  the  vertical  and  horizontal  position  tolerances  for  a 
particular  angle  of  incidence.  These  tolerances  can  then  be  compared  with  the 
vertical  and  horizontal  position  resolution  of  an  aircraft  tracking  system. 
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Inspecting  F  igures  1-14  and  1-15,  the  following  observations  may  be 

made: 

o  Near  overhead  (zero  angle  of  incidence)  position  tolerances  are 
large,  both  for  vertical  and  horizontal  deviations. 

o  For  vertical  deviations,  position  tolerance  rapidly  decreases  with 
increasing  angle  of  incidence,  and  tends  toward  0  at  grazing 
incidence. 

o  For  horizontal  deviations,  the  angles  of  incidence  most  sensitive  to 
reflection  correction  are  around  35  +  (10  to  20)  degrees,  the 
interval  increasing  with  decreasing  altitude. 

As  an  example,  consider  the  Aircraft  Radar  Tracking  System  (ARTS  III) 
used  by  FAA  at  many  busy  airports.  Its  nominal  vertical  resolution  is  +  50  feet  =  15 
meters,  its  nominal  horizontal  resolution  is  +  0.005  nautical  miles  =  +  10  meters. 
Choosing  e  =  0. 1  dB  we  can  see  that  for  an  altitude  of  100  m,  angles  of  incidence  of 
+  40°  will  have  less  than  the  specified  error.  That  angle  interval  increases  to  +  65° 
for  an  altitude  of  500  m.  It  appears,  therefore,  that,  at  least  for  overhead  flyover, 
the  ARTS  system's  resolution  should  be  sufficient  for  defining  an  aircraft's  position 
for  applying  the  reflection  correction.  FHowever,  there  exists  other  problems 
concerning  application  of  a  reflection  correction  which  are  further  discussed  in 
Section  l. 3.2.3. 


1.3.4  Frequency  Cutoff  Technique 

This  technique  qttempts  to  compensote  for  the  occurrence  of  pseudo¬ 
tones  by  restricting  the  range  of  frequencies  within  which  tone  corrections  may  be 
calculated,  while  PNL  itself  is  calculated  over  the  entire  frequency  range. 
Currently,  FAR  Part  36  allows  cutting  off  any  tone  corrections  below  800  Hz  if  it 
can  be  demonstrated  that  there  are  no  real  tones  below  that  cutoff  which  might 
cause  a  greater  tone  correction  than  any  tones  above  the  cutoff. 

In  order  to  investigate  the  effect  of  various  cutoff  frequencies,  many 
pseudotone  contaminated  spectral  time  histories  from  the  aircraft  flyovers 
recorded  at  LAX  were  subjected  to  an  EPNL  analysis  allowing  the  cutoff  to  take  on 
the  following  values:  50  Hz  (no  cutoff),  400  Hz,  800  Hz,  1600  Hz. 

Table  1-3  shows  the  results  listing  the  differences  between  EPNL's  with 
and  without  tone  correction  cutoff.  Nonzero  differences  indicate  that  for  at  least 
a  portion  of  the  spectrum  within  the  10-dB-down  interval,  there  exist  tones  in  the 
lower  part  of  the  spectrum,  may  they  be  pseudotones  or  real  tones.  Each  individual 
flyover  time  history  is  identified  by  what  is  called  a  dataset  sentinel  which  refers 
to  corresponding  spectral  time  history  plots  in  Appendix  A.  Note  that  the  datasets 
there  are  arranged  in  a  sequence  given  by  Table  1-7  to  be  found  in  Section  1.4.3. 3. 
Some  approach  flights  occurred  displaced  laterally  by  about  230  m  (750  ft)  from 
the  overhead  position.  These  are  identified  as  "sideline"  in  Table  1-3.  These 
fliqhts  are  not  representative  of  the  FAR  Part  36  approach  measurement  con¬ 
ditions,  but  they  are  included  as  a  matter  of  interest. 

With  very  few  exceptions,  the  EPNL  difference  magnitudes  are  less  than 
O.b  dB,  Since  this  is  within  expected  measurement  errors,  one  would  be  tempted  to 
dispense  with  the  subject  and  conclude  that  an  application  of  frequency  cutoff  for 
tone  correction  has  a  negligible  influence.  However,  during  certification  of  large 
passenger  aircraft,  minute  differences  can  decide  passing  or  failing  of  FAR  Part  36 
requirements.  Differences  in  EPNL  values  of  the  order  of  0. 1  dB  are  therefore 
important  to  the  aircraft  noise  certification  process. 

Analyzing  the  results  by  aircraft  type,  it  may  be  seen  that,  for  the  B707 
(representing  the  four-engine  narrow  body  low  engine  bypass  ratio  aircraft  cate¬ 
gory),  EPNL  difference  magnitudes  are  mostly  zero  or  very  small.  This  indicates 
the  presence  of  strong  tones  above  any  tone  correction  cutoff,  even  for  the  takeoff 
condition  where  one  might  expect  jet  noise  to  dominate  and  mask  tones.  However, 
even  if  such  masking  did  occur,  the  tone  correction  as  presently  prescribed  by  FAR 
Part  36  would  still  be  added  at  full  value  (this  appears  to  be  a  weakness  of  the  tone 
correction  procedure). 
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Table  1-3 


Differences  Between  EPNL's  Calculated  with  Various  Tone  Correction 
Lower  Cutoff  Frequencies 
Ground  Surface  Hard  Unless  Otherwise  Noted 
EPNL  Subscript  Indicates  Lower  Cutoff  Frequency  (Hz)  for  Tone  Correction 


Microphone 

Location 

Aircraft 

Dataset 

Sentinel 

EPNLjq  -  EPNL^jq 

dB 

EPNLsq  -  EPNL^ 

dB 

EPNLso-EPNL^oo 

dB 

Approach 

1.2  m 

B727 

(Soft  Surfoce) 

19010 

0.04 

0.04 

0.04 

B727 

21010 

0.00 

0.00 

0.00 

B727 

(Soft  Surface) 

42010 

0.06 

0.06 

0.06 

B727 

50010 

0.08 

0.10 

0.16 

B707 

(Soft  Surface) 

41010 

0.00 

0.00 

0.00 

"Sideline" 

B707 

(Soft  Surface) 

43010 

0.00 

0.00 

0.00 

B707 

49010 

0.00 

0.00 

0.00 

"Sideline" 

B707 

51010 

0.00 

0.00 

0.00 

"Sideline" 

DC-10 

14010 

0.00 

0.28 

0.31 

"Sideline" 

DC- 10 

(Soft  Surfoce) 

17010 

0.00 

0.34 

0.38 

DC- 10 

48010 

(0.24)* 

(0.83) 

(1.06) 

Sideline, 

1.2  m 

B727 

29010 

0.02 

0.24 

0.26 

B727 

64020 

0.07 

0.24 

0.43 

B707 

65020 

0.00 

0.00 

0.00 

B707 

26010 

0.00 

0.00 

0.00 

B707 

27010 

0.02 

0.06 

0.11 

DC- 10 

66010 

(0.02) 

(0.03) 

(0.04) 

DC- 10 

67010 

0.07 

0.11 

0.11 

Table  1-3  (Continued) 


Microphone 

Dataset 

EPNLS0  -  EPNL^q 

EPNL^  -  EPNLgQQ 

EPNL50  -  EPNL  ,600 

Location 

Aircraft 

Sentinel 

dB 

dB 

dB 

Sideline, 

B727 

28010 

0.08 

0.08 

0.13 

10  m 

B727 

60020 

0.36 

0.37 

0.38 

B707 

24010 

0.00 

0.00 

0.20 

B707 

25010 

0.08 

0.08 

0.29 

B707 

61020 

0.00 

0.00 

0.00 

DC-10 

62020 

(0.08) 

(0.08) 

(0.08) 

DC- 10 

63020 

0.09 

0.16 

0.16 

Takeoff, 

B727 

36010 

0.20 

0.54 

0.54 

1.2  m 

B727 

92020 

0.01 

0.20 

0.25 

B707 

94020 

0.00 

0.00 

0.00 

B707 

9S020 

0.00 

0.00 

0.00 

DC-10 

96020 

0.06 

0.12 

0.12 

DC- 10 

97020 

0.13 

0.23 

0.24 

8747 

91020 

0.00 

0.00 

0.00 

B747 

93020 

0.00 

0.00 

0.00 

*  Numbers  in  parentheses  are  derived  from  EPNL  values  calculated  from  incomplete  spectral 
time  histories  (i.e,  at  least  one  of  the  10-dB-down  points  of  the  PNLT  time  history  was  not 
reached). 
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The  B727  represents  the  JT8D-powered  aircraft.  The  trend  for  the 
EPNL  differences  here  is:  very  small  for  approach,  larger  for  sideline  and  takeoff. 
This  is  more  of  an  expected  trend:  on  approach  the  high-pitch  fan  tones  dominate 
the  tone  correction;  on  takeoff,  pseudotones  may  cause  greater  tone  corrections 
than  real  tones. 

The  DC- 10  and  B747  represent  the  wide  body  high  bypass  ratio  aircraft. 
The  B747  takeoffs  are  obviously  tone  correction  dominated  at  higher  than  the 
cutoff  frequencies.  For  the  DC- 10,  EPNL  differences  are  generally  larger  than  for 
other  aircraft  types.  This  could  either  be  because  there  exist  real  tones  in  the 
lower  part  of  the  spectrum,  or  due  to  pseudotones.  Knowing  the  noise  emission 
characteristics  of  a  high  bypass  ratio  engine,  the  likely  situation  is  a  combination 
of  both  real  and  pseudotones.  However,  inspection  of  the  spectral  time  history 
plots  of  Appendix  A  indicates  that  pseudotones  are  mostly  responsible  for  low 
frequency  tone  corrections. 

It  is  difficult  to  derive  general  recommendations  from  the  above  results. 
For  overhead  flyovers  (approach,  takeoff)  the  method  of  restricting  the  frequency 
range  above  of  where  pseudotones  occur  does  work  quite  well  providing  there  are 
no  real  tones  below  that  range.  For  sideline,  the  angle  of  incidence  is  usually  so 
large  as  to  push  the  pseudotones  up  into  the  vicinity  of  real  fan,  turbine,  or 
compressor  tones  if  the  1.2  m  microphone  height  is  used.  In  this  case,  a  frequency 
cutoff  technique  cannot  be  used.  For  the  10  m  microphone  and  sideline, 
pseudotones  occur  in  about  the  same  frequency  range  as  for  the  1.2  m  microphone 
for  overhead  flyovers,  so  that  this  technique  may  again  be  applied.  However,  there 
are  other  reasons  that  make  this  frequency  cutoff  method  unattractive: 

o  The  effects  of  pseudotones  are  not  corrected  for,  but  simply 
"evaded"  as  far  as  the  tone  correction  is  considered.  PNL  still  is 
calculated  in  the  presence  of  the  pseudotones. 

o  A  narrow  band  analysis  (not  done  in  this  study)  must  be  performed 
to  prove  the  absence  of  real  tones  below  the  cutoff  frequency. 

o  The  method  cannot  be  applied  generally  for  all  aircraft  types  such 
as  propeller  and  rotor  aircraft  which  usually  exhibit  strong  tones  in 
the  lower  frequency  regions. 

Section  1.5  attempts  to  strike  a  balance  between  the  various  pseudotone 
correction  procedures  considered  in  this  study,  and  incorporates  the  results  of  this 
section. 
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1.4  Instrumentation  Correction  Techniques 

We  will  define  instrumentation  correction  techniques  for  pseudotones  as 
those  procedures  which  are  applied  either  by  physical  arrangement  of  the  measure¬ 
ment  setup,  or  are  applied  to  the  instantaneous  acoustic  signal  before  it  enters  a 
frequency  decomposition  device.  These  solutions  are  attractive  because  no 
analytical  corrections  will  be  necessary  when  the  data  is  reduced.  The  raw  data 
will  reflect  only  noise  due  to  the  source  (aircraft)  and  therefore  should  be  less 
dependent  on  the  measurement  site  than  data  contaminated  by  reflections. 

As  discussed  in  Section  1. 1,  the  pseudotones  exist  due  to  an  interaction 
of  the  incident  noise  with  the  ground  plane.  The  actual  pseudotones  which  occur 
are  a  function  of  the  impedance  of  the  ground  plane,  the  position  of  the  aircraft  in 
relation  to  the  microphone,  and  the  height  of  the  microphone  above  the  ground 
plane.  Alternate  measurement  techniques  could  be  designed  to  alter  one  or  more 
of  these  parameters  in  order  to  reduce  or  eliminate  the  effect  of  pseudotones. 
Possible  techniques  are: 

o  Surface  or  flush-mounted  microphone 

o  Elevated  microphone 

o  Simulated  "anechoic"  ground  plane 

o  Highly  directional  microphone 

o  Moving  microphone  array 

o  Real  time  electronic  cancellation  system 

o  Real  time  impedance  cancellation  system 

1.4.1  Overview  of  Methods  (See  Table  1-4) 

a.  Surface  or  Flush-Mounted  Microphone 

The  microphone  is  either  mounted  directly  into  the  ground  surface  so 
that  its  diaphragm  is  in  the  same  place  as  the  surface,  or  it  is  laid  on  the  surface 
with  the  diaphragm  perpendicular  to  the  ground  surface,  or  it  is  pointed  down 
towards  the  ground  surface  and  fixed  at  a  very  small  distance  from  it  (usually  1/2 
inch).2’  8 
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This  technique  nearly  eliminates  the  reflected  wave  problem  (see  Section 
1.4.3. 1  for  a  further  discussion  on  expected  performance  of  a  ground  microphone). 
Pressure  doubling  occurs  at  the  ground  surface  (assuming  it  is  rigid)  which 
increases  the  signal  by  6  dB  relative  to  free  field.  A  series  of  spectra  obtained 
with  this  technique  is  shown  in  Figure  1-16  which  should  be  compared  with  the 
earlier  Figure  1-3  showing  spectra  measured  simultaneously  with  a  1.2  m  micro¬ 
phone.  The  absence  of  pseudotones  for  the  flush-mounted  microphone  is  striking. 

This  is  a  very  attractive  technique  because  it  is  relatively  simple  to  set 
up  and  the  microphone  is  easy  to  get  at  for  calibration.  In  addition,  data  should  be 
compatible  between  measurement  sites.  However,  the  effects  of  surface  smooth¬ 
ness  and  size  as  well  as  the  often  strong  refractive  effects  of  wind  and 
temperature  gradients  very  near  the  surface  are  not  well  understood. 

b.  Elevated  Microphone 

The  pseudotones  which  are  observed  in  aircraft  flyover  noise  signatures 

are  a  function  of  the  height  of  the  microphone  above  the  reflecting  surface.  If  the 

height  of  the  microphone  is  increased  to  10  meters  or  so,  pseudotones  are  pushed 

down  in  frequency  to  a  range  we  are  not  interested  in  (except  for  sideline;  see 

Figure  1-25  in  Section  1.4. 3.1).  This  technique  would  allow  measurements  to  be 

taken  over  almost  any  ground  surface.  The  main  drawback  of  this  technique  is  that 

access  to  the  microphone  is  awkward.  Calibration  and  setup  are  less  convenient, 

although  manageable.  In  addition,  a  10  meter  or  higher  structure  could  be 

undesirable  near  an  airport.  Nevertheless,  one  major  aircraft  engine  manufacturer 

9  23  27 

strongly  advocates  this  technique.  ’  ' 

c.  Anechoic  Ground 

This  technique  is  now  being  used  for  some  static  aircraft  engine  testing 
and  it  involves  applying  absorptive  material  or  an  anechoic  section  consisting  of 
sound-absorbing  wedges  on  the  ground  surface  where  the  reflections  would  occur.  ^ 
This  method  eliminates  the  reflected  wave  and  hence  the  pseudotones  and  also 
maintains  the  1.2  meter  microphone  height  which  is  convenient  for  access  and 
above  the  strongest  flow  and  thermal  gradients  near  the  ground.  However,  for 
flyover  measurements,  the  amount  of  absorptive  material  required  may  be  prohi¬ 
bitive,  or  the  anechoic  section  would  have  to  move  with  the  aircraft  position. 


igure  1-16.  Flyover  Spectra  at  1/2  Second  Intervals,  Ground  Level  Microphone 
(from  Reference  2) 


d.  Highly  Directional  Microphone 


If  a  highly  directional  microphone  is  pointed  at  the  aircraft  during  the 
flyover,  all  but  the  direct  sound  will  be  virtually  eliminated.  This  technique  has  an 
additional  benefit  in  that  acoustic  ambient  noise  will  be  reduced.  A  method  of 
tracking  the  aircraft  with  sound  travel  time  delay  would  be  necessary,  and  since 
the  frequency  response  of  directional  microphones  is  usually  nonuniform,  a 
frequency  correction  would  have  to  be  applied  during  data  reduction.  Although 
current  technology  is  available  for  implementing  this  approach,  including  the  use  of 
electronically  steerable  arrays,  it  is  not  a  practical  approach  for  general  applica¬ 
tion  at  this  time. 

e.  Moving  Microphone  Array 

Microphones  at  differing  heights  which  are  multiplexed  to  obtain  an 
average  sound  level  will  produce  results  in  which  the  pseudotones  are  averaged  and 
hence  greatly  reduced.  Moving  the  microphones  will  reduce  the  probability  that  a 
pseudotone  at  a  given  frequency  will  continue  to  occur.  This  method  would  reduce 
the  influence  of  ground  reflections  and  make  data  more  repeatable  between 
measurement  sites.  There  are  several  development  problems  associated  with  this 
method,  one  of  which  is  designing  a  mechanical  system  which  will  move  the 
microphones  without  producing  mechanical  noise  or  developing  an  electronic 
equivalent  noise.  Neglecting  this  substantial  problem,  considerable  effort  would 
also  be  required  to  standardize  the  dimensions  and  effective  motion  of  the  array. 

f.  Real  Time  Electronic  Cancellation  System 

This  method  would  include  a  microprocessor  which  would  be  fed  data  on 
the  distance  and  position  of  the  aircraft  in  order  to  calculate  instantaneous 
reflection  corrections  and  thereby  suppress  pseudotones  at  the  input.  This  highly 
sophisticated  system  would  be  difficult  to  develop.  However,  once  developed,  it 
would  produce  very  repeatable  data.  The  system  could  be  conveniently  calibrated 
by  feeding  it  a  simulated,  reflected  wave  only  and  checking  +hat  the  output  is  zero. 
Since  it  would  use  reflection  theory,  this  technique  suffers  from  some  of  the 
drawbacks  listed  for  the  analytical  reflection  correction  technique;  in  particular, 
surface  impedance  should  be  known  as  accurately  as  possible.  Furthermore,  the 
general  lack  of  success  in  this  program  of  applying  an  analytically-derived 
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correction  in  the  data  processing  phase  argues  against  use  of  this  similar  but  real 
time  approach. 


1.4.2  Selection  of  Instrumentation  Correction  Techniques 

The  same  criteria  as  listed  at  the  beginning  of  Section  1.3.1  are  used  for 
selection  of  an  instrumentation  technique  to  eliminate  or  minimize  pseudotones. 

The  technical  feasibility  of  the  surface  microphone  has  been  demon¬ 
strated.  The  accuracy  of  the  method  needs  to  be  investigated.  It  is  cheaply 
implemented,  easily  enforceable  because  appropriate  specifications  can  be  written 
straightforwardly,  and  it  could  be  implemented  on  relatively  short  notice. 

The  technical  feasibility  of  the  elevated  microphone  technique  has  been 
demonstrated;  its  accuracy  is  controversial.  It  is  not  cheaply  implemented  because 
of  the  tower  requirement  and  the  awkward  access  to  the  microphone.  The 
technique  would  be  easily  enforceable  and  could  be  implemented  on  relatively  short 
notice. 

The  anechoic  ground  technique  is  technically  feasible  and  very  accurate 
if  practically  all  reflections  can  be  eliminated.  It  would  be  quite  expensive  as  an 
outdoor  rig  would  be  required  which  effectively  absorbs  sound  over  most  of  the 
audible  range  and  which  would  have  to  be  very  large  (or  moving  with  aircraft 
position)  to  cover  all  possible  reflection  points.  Enforcement  of  this  technique 
could  turn  into  an  administrative  nightmare  considering  that  the  absorptive 
qualities  of  the  surface  would  need  to  be  demonstrated.  The  time  frame  of 
implementation  would  probably  be  several  years  including  the  development  and 
construction  of  a  suitable  surface. 

The  use  of  a  highly  directional  microphone  would  be  technically  feasible 
and  probably  very  accurate.  The  technique  would  be  expensive  as  the  microphone 
would  have  to  be  pointed  in  the  direction  of  sound  arrival  which  is  different  from 
the  optical  aircraft  position;  some  kind  of  computerized  pointing  system  is 
probably  required.  Because  of  this  complicated  system,  enforcement  would  be 
difficult,  and  it  would  take  a  long  time  to  implement. 

The  real  time  electronic  cancellation  system  may  be  technically  feasible 
with  a  state-of-the-art  digital  processor  modeling  the  time-domain  transfer 


function  of  the  reflection  process.  It  would  be  very  expensive,  probably  difficult  to 
enforce,  and  would  take  a  long  time  to  implement. 

The  above  discussion  clearly  shows  that  there  are  two  rather  straight¬ 
forward  instrumentation  correction  techniques  for  pseudotones  which  show  good 
promise  of  being  successfully  applied  in  routine  aircraft  noise  certification 
measurements.  These  techniques  are  the  flush-mounted  and  elevated  microphones. 
These  techniques  are  investigated  in  more  detail  in  the  subsequent  sections  of  this 
report. 

1.4.3  Surface  and  Elevated  Microphones 

Smith  (Reference  23)  and  McKaig  (Reference  2)  have  given  excellent 
reviews  of  the  uses  and  effects  of  surface  and  elevated  microphones,  both  also 
reporting  measured  data.  Here,  their  results,  conclusions  and  recommendations  are 
combined  with  theoretical  considerations  and  with  experimental  evidence  collected 
under  this  study. 

1.4.3. 1  Theoretical  Considerations 

It  is  instructive  to  study  first  some  theoretical  results  obtained  by 
applying  the  reflection  algorithm  given  in  Section  1.3. 2. 1  to  the  surface  and 
elevated  microphones.  A  simple  way  to  mount  a  surface  microphone  is  to  invert  it 
and  position  the  diaphragm  at  a  very  small  distance,  say,  1/2  inch,  from  the  ground 
surface.  Under  the  assumption  of  an  infinitely  large  and  plane  surface  of  large 
impedance  (flow  resistivity  10^  cgs  rayls/cm),  and  a  point  microphone,  Figure  1-17 
shows  what  correction  would  need  to  be  added  to  the  measured  signal  in  order  to 
obtain  free  field  conditions.  It  is  seen  that,  particularly  near  the  overhead  position 
(small  angle  of  incidence),  that  correction  is  not  even  near  the  expected  -6  dB  for 
high  frequencies  (4  kHz  and  up)  often  strongly  contributing  to  PNL. 

This  problem  could  be  circumvented  by  mounting  the  microphone  truly 
flush  with  the  ground  surface.  This  usually  necessitates  mounting  the  microphone 
through  a  board  or  panel.  The  minimum  size  requirements  for  such  a  panel  can  be 
estimated  from  diffraction  theory.  Figure  1-18  is  taken  from  Reference  24.  It 
indicates  that  the  product  of  wave  number  K  and  disc  radius  a  should  be  at  least 
10  in  order  to  acoustically  simulate  an  infinite  plane.  This  translates  into  the 
following  requirement: 
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Figure  1-17.  Plot  of  Spectral  Time  History  of  Reflection  Correction  for  a  Simulated 
Takeoff  Condition:  Climb  Angle  66,  Altitude  at  Overhead  244  m, 
Speed  80  m/sec,  Microphone  Height  1/2  inch. 
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disc  radius  should  be  greater  than  (10  times  speed  of  sound)  /  (2TTf) 

where  f  is  frequency.  The  lowest  frequencies  of  interest  are  around  50  Hz.  With 
343  m/sec  for  the  speed  of  sound,  the  minimum  disc  diameter  becomes  about  22  m 
(about  70  feet),  a  very  large  disc  indeed.  Even  if  requirements  could  be  relaxed  by 
a  factor  of  2  or  3  (by  strategically  placing  the  microphone  on  the  panel,  giving  the 
panel  a  noncircular  shape,  raising  the  lowest  frequency  of  interest),  the  panel  size 
would  still  be  quite  large.  At  permanent  stationary  test  installations  this  may 
quite  easily  be  accomplished  by  imbedding  a  small  panel  flush  info  a  larger  hard 
surface  (asphalt  or  concrete).  A  portable  system  seems  indeed  difficult  to  realize. 

Turning  towards  the  elevated  microphone,  let  us  compare  Figures  1-8, 
1-19,  and  1-20,  which  demonstrate  (for  a  soft  surface)  how  the  interference 
pattern  moves  to  lower  frequencies  as  the  microphone  height  is  increased  from  1.2 
m  to  5  m  and  10  m.  The  pattern  in  the  vicinity  of  the  overhead  position  is  of 
course  most  important  since  noise  levels  are  greatest  there.  Whereas  the  5  m 
location  still  shows  significant  corrections  above  50  Hz  (8-th  and  9-th  spectra  from 
left  of  Figure  1-19  are  closest  to  overhead),  the  effect  at  10  m  can  be  termed 
almost  negligible  at  overhead  and  is  conditionally  acceptable  further  away  from 
overhead;  the  word  "conditionally"  is  inserted  because  the  reflection  effects  may 
be  unacceptable  in  the  case  of  a  very  directional  source  where  the  strongest 
flyover  noise  level  occurs  away  from  overhead,  possibly  in  a  region  where  even  the 
pseudotones  for  a  10  m  microphone  height  may  not  have  a  significant  influence  on 
EPNL. 

Figure  1-21  portrays  the  same  case  as  Figure  1-20  except  that  the 
acoustic  surface  impedance  is  very  large  simulating  a  rigid  surface.  At  low 
frequencies  there  is  hardly  any  difference  between  soft  and  hard  ground,  i.e.,  soft 
grass  acts  like  a  hard  surface.  At  higher  frequencies  the  soft  surface  absorbs  much 
or  most  of  the  sound  (approximately  0  dB  correction)  whereas  the  hard  surface 
reflects  it  (about  -3  dB  correction  to  free  field). 

The  preceding  discussion  applied  to  overhead  flyovers.  Both  the  ground 
microphone  and  elevated  microphone  are  reasonably  free  of  pseudotones.  The 
elevated  microphone  also  avoids  possible  panel  diffraction  problems  which  may 
occur  with  a  casual  ground  microphone  installation.  At  sideline  conditions,  the 
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theoretically  predicted  interference  patterns  look  very  different.  Figure  1-22 
shows  the  reflection  correction  for  a  simulated  takeoff:  climb  angle  is  9°,  at  the 
sideline  measurement  point,  the  elevation  angle  is  20°  and  the  aircraft  altitude  and 
speed  are  200  m  and  80  m/sec,  respectively.  Ground  was  assumed  soft  with  a  flow 
resistivity  of  300  cgs  rayls/cm;  the  associated  impedances/admittances  as  a 
function  of  frequency  are  shown  in  Table  1-5.  Microphone  height  is  1.2  m.  Near 
grazing  incidence  (90°  angle  of  incidence,  0°  elevation  angle),  the  correction  shows 
a  broad  hump  centered  around  400  Hz.  This  hump  begins  to  clearly  develop  at 
elevation  angles  of  about  5°  and  below.  Such  behavior  can  be  observed  in  measured 
spectra  which  show  a  marked  deficiency  in  the  same  frequency  region  for  sound 
which  has  traveled  horizontally  over  acoustically  soft  ground.  The  latter  part  of 
the  spectral  time  history  of  Figure  1-2  shows  the  pseudotone  interference  pattern 
as  has  been  observed  before.  The  pattern  in  the  early  part,  however,  moves  up  to 
higher  frequencies  as  should  be  expected  from  the  reduced  path  length  difference 
between  direct  and  reflected  waves  (note  that  this  implies  that  the  analytical 
pseudotone  correction  technique  consisting  of  disregarding  tone  corrections  to  PNL 
below  a  certain  frequency  may  fail  on  sideline). 

The  input  to  Figure  1-23  differs  from  that  of  Figure  1-22  merely  by  the 
nature  of  the  ground  which  is  assumed  hard  with  a  flow  resistivity  of  10'*  cgs 
rayls/cm.  The  resulting  impedances/admittances  are  shown  in  Table  1-6.  Note 
that  the  grazing  incidence  hump  develops  at  much  higher  frequencies.  This  case  is, 
however,  hardly  of  practical  significance  as  sideline  propagation  will  most  fre¬ 
quently  occur  over  acoustically  soft  ground. 

The  spectral  time  history  plot  for  the  reflection  correction  at  sideline 
conditions  for  the  surface  microphone  (assumed  to  be  truly  flush  mounted  in  an 
infinite  acoustically  soft  surface)  does  not,  of  course,  show  any  pseudotone  pattern, 
but  it  also  is  not,  by  any  means,  constant.  A  correction  to  free  field  would  have  to 
account  for  ground  absorption  which  is  what  Figure  1-2 4  displays.  If  one  were  to 
place  a  large  hard  disc  (10  m  diameter  or  so)  around  the  microphone,  the 
theoretical  handling  of  the  case  of  near  grazing  incidence  would  be  made  very 
difficult  or  impossible  since  sound  waves  traveling  horizontally  (or  almost  so)  will 
encounter  a  surface  impedance  discontinuity.  This  analysis  therefore  suggests 
that,  if  a  surface  microphone  is  at  all  used  for  sideline  measurements,  it  be  placed 
right  on  the  soft  surface.  An  attempt  could  then  be  made  to  correct  for  ground 
absorption. 
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Table  1-5 


Impedance/Admittance  for  a  Surface  Flow  Resistivity  of  300  cgs  rayls/cm 
Imaginary  Part  of  Impedance  Is  Negative  (See  Footnote  on  Page  1-22) 
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SPECTRAL  PL0T  F0R  CONSECUTIVE  TIME  INTERVALS 


Table  1-6 


Impedance/Admittance  for  a  Surface  Flow  Resistivity  of  10  cgs  rayls/cm 
Imaginary  Part  of  Impedance  Is  Negative  (See  Footnote  on  Page  1-22) 
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Figure  1-24.  Ground  Effect  Correction  for  Sideline,  Soft  Ground,  Ground 
Microphone.  Elevation  Angles  Some  as  on  Figure  1-22. 
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Figure  1-25  shows  the  calculated  ground  interference  pattern  for  side¬ 
line  with  an  elevated  microphone  (10  m  above  ground).  The  figure  shows  that 
pseudotone  phenomena  occur  near  the  closest  approach  portion  of  the  flyby  (middle 
of  plot).  In  particular,  the  inverted  peak  at  100  Hz  in  the  six  lower-most  spectra 
could  give  rise  to  a  spurious  tone  correction.  At  very  low  elevation  angles  (top 
portion  of  plot)  the  pattern  shifts  to  higher  frequencies  and  the  above  mentioned 
"hump"  appears  at  lower  frequencies. 

The  above  paragraphs  on  sideline  measurements  indicate  that  it  is  very 
difficult  to  obtain  free  field  spectra  from  sideline  measurements  considering  only 
the  presence  of  an  idealized  ground.  This  is  compounded  by  several  other  effects 
which  introduce  largely  unpredictable  variations: 

o  Terrain  irregularities  and  obstacles  (diffraction) 

o  Temperature  and  wind  gradients  (refraction) 

o  Atmospheric  turbulence  which  is  strongest  near  the  ground  (scat¬ 
tering) 

Sideline  measurements  are  much  more  susceptible  to  the  above  influences  than 
approach  or  overhead  takeoff  measurements  since  the  propagation  path  is  much 
closer  to  the  ground  (much  smaller  elevation  angles  during  high-noise  portion  of 
flyby)  for  the  sideline  condition  than  for  the  other  two  measurement  locations. 

1. 4.3.2  Short  Literature  Review 

McKaig  (Reference  2)  states  that  "...it  appears  hopeless,  or  at  least 
extremely  risky,  to  try  to  compensate  analytically  for  the  presence  of  ground 
reflection."  He  reports  on  the  use  of  surface  microphones  1/2  inch  from  the 
surface  saying  that  at  about  14  kHz  there  will  be  at  most  a  0.2  dB  error  from  the 
6  dB  perfect  reinforcement  case  with  the  difference  being  less  at  lower  fre¬ 
quencies.  This  is  not  at  all  in  agreement  with  the  results  shown  in  Figure  1-17 
even  considering  that  the  finite  microphone  size  will  reduce  the  theoretical 
corrections  of  that  figure.  McKaig  also  found  that  near-grazing  incidence  angles 
did  not  produce  any  loss  of  high  frequency  signal  strength  so  long  as  the  angle  of 
incidence  was  less  than  87°  (elevation  angle  greater  than  3°).  He  recognizes  that 
the  ground  microphone  technique  may  not  be  at  all  applicable  to  sideline 
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measurements.  Use  of  flush-mounted  microphones  in  acoustic  scale  model  tests 
does  not  permit  the  customarily  casual  full-scale  mounting  techniques;  good  results 
were  obtained  only  with  extreme  care,  figures  1-3  and  1-16  presented  earlier 
were  taken  from  McKaig's  paper. 

In  Reference  23,  Smith  makes  an  excellent  case  against  the  1.2  m 
microphone  height,  and  in  favor  of  the  elevated  microphone  as  opposed  to  the 
ground  microphone.  Some  pertinent  quotations  from  Smith's  paper: 

o  "1.2  and  10  meter  locations  produce  virtually  the  same  absolute 
levels  over  varied  surfaces,  except  at  low  angles  to  the  flight  line 
where  the  1.2  meter  results  are  totally  distorted.  A  flush  system 
shows  amplification  relative  to  the  other  systems,  3  dB  to  5  dB. 

o  Surface  conditions  have  a  similar  effect  on  the  1.2  and  10  meter 
systems.  The  impedance  change  from  board  to  grass  with  the  flush 
system  appears  to  produce  some  spectral  change  relative  to  an 
infinite  hard  surface,  which  would  be  difficult  to  allow  for  in 
practical  tests. 

o  At  1.2  meters,  ground  interference  effects  produce  extreme  spec¬ 
tral  distortion,  causing  loss  of  several  1/3  octave  bands  at  shallow 
angles  of  incidence.  10  meter  and  flush  systems  minimize  or 
eliminate  interference  effects. 

o  Recognizable  source  characteristics  are  well  displayed  by  both  a 
flush  system  in  a  hard  surface  and  by  10  meter  systems  over  both 
hard  and  soft  surfaces.  Source  character  is  unacceptably  distorted 
by  a  1.2  meter  system. 

o  A  system  mounted  10  meters  above  ground  produces  less  scatter 
due  to  ground  surface  conditions,  less  destructive  impact  from 
ground  reflection  and  a  more  faithful  portrayal  of  source  char¬ 
acter.  In  a  certification  context  this  virtually  eliminates  "pseudo 
tone"  corrections,  avoids  the  effect  of  localized  temperature  and 
humidity  gradients  and  does  not  introduce  any  increased  back¬ 
ground/wind  noise.  Such  a  system  also  minimizes  the  problem  of 
1/3  octave  "drop  out"  at  the  fringes  of  the  time  history. 
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o  A  change  to  a  10  meter  location  would  put  the  measuring  systems 
at  the  same  position  that  is  currently  used  to  define  ground  level 
atmospheric  conditions." 

A  surprising  statement  is  that  the  use  of  the  elevated  microphone  does 
not  introduce  any  increased  background/wind  noise  since  the  opposite  might  be 
expected  as  wind  speed  generally  increases  with  height  above  ground.  The 
explanation  probably  is  that  l(J  m  is  usually  just  above  the  atmospheric  viscous 
boundary  layer  with  its  turbulence  adversely  affecting  lower  microphones  whereas 
the  10  m  microphone  is  more  likely  to  be  exposed  to  a  steady  stream  whose  speea 
may  be  less  than  the  peak  velocities  in  the  turbulence  below. 

Smith  has  also  drafted  an  Appendix  to  a  proposed  SAE  AIR  (Refer¬ 
ence  27)which  is  based  on  the  same  measurements  as  Reference  23,  but  presents 
some  additional  relevant  details  and  observations: 

o  Solar  heating  and  small  scale  turbulence  above  a  surface  mounted 
microphone  cause  high  frequency  loss;  the  turbulence  also  gener¬ 
ates  "comparatively  high  levels  of  ambient  noise." 

o  Diffraction  ripples  induced  by  the  edge  of  an  8  feet  diameter  board 
used  for  mounting  the  microphone  are  generally  within  3  to  4  dB. 

o  At  the  10  m  microphone  location,  there  are  some  small  "augmenta- 
tional/cancellation"  effects,  i.e.,  recognizable  reflection  inter¬ 
ference  patterns,  in  the  very  low  frequencies.  This  corroborates 
the  theoretical  prediction  of  Figure  1-20. 

Smith  recommends  that  microphones  mounted  10  m  above  a  hard  ground 
be  used.  The  surface  should  be  painted  a  light  color  to  avoid  solar  heating  effects. 
If  the  measurements  are  taken  over  acoustically  soft  terrain,  he  proposes  a 
correction  curve  (shown  in  Figure  1-26)  replacing  the  constant  -3  dB  correction  to 
be  applied  for  hard  surfaces. 

Reference  25  demonstrated  that  at  a  height  of  137  m  (450  feet)  above 
the  ground,  aircraft  noise  spectra  are  completely  free  of  ground  reflections. 
Reference  26  found  that  a  reflection  signal  comparable  to  the  diret  ‘  signal  can  be 
found  at  heights  of  up  to  60  m  (200  feet)  above  the  ground  for  frequencies  of  less 
than  1000  Hz. 
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1.4.3. 3  Evaluation  of  Data  Measured  at  LAX 


Appendix  A  contains  Spectral  Time  History  (STH)  plots  (spectra  at 
successive  time  intervals)  and  for  various  microphone  locations  and  ground  surfaces 
employed  during  the  measurement  of  over  20  flights  from  Los  Angeles  Inter¬ 
national  Airport.  Table  1-7  summarizes  the  observations  which  may  be  made  on 
inspection  of  these  STHs.  The  table  also  lists  numerical  values  for  each  flight. 
The  "waviness"  is  defined  for  a  single  spectrum  as  follows: 

32 

w=  I  <LW  I  -  2Li  *  Li  - 1>2 
i  =  18 

where  Lj  is  the  level  of  the  i-th  third-octave  band,  with  i  indicating  standard  third- 
octave  band  numbers.  W  is  a  measure  of  how  much  the  spectrum  fluctuates  in  the 
frequency  domain.  Mathematically,  W  is  the  sum  of  squares  of  numerical 
approximations  to  the  second  derivative  of  the  spectrum  with  respect  to  the 
logarithm  of  frequency.  This  calculation  is  limited  to  the  lower  portion  of  each 
spectrum  (below  2  kHz)  where  pseudotones  commonly  occur  since  W  is  used  as  a 
numerical  indicator  of  the  power  of  correction  methods  to  remove  pseudotones. 
Two  waviness  values  are  listed  in  Table  1-7,  one  for  the  spectrum  at  PNLTm(JX, 
and  the  other  is  the  mean  +  one  standard  deviation  over  the  10-dB-down  interval. 

The  noise  levels  listed  in  Table  1-7  were  computed  in  the  standard 
manner.  The  integral  of  the  A-weighted  time  history  is  the  sound  exposure  level 
L(~  for  which  the  integration  time  interval  was  taken  as  the  10-dB-down  interval. 
The  latter  may  be  different  from  the  one  used  for  EPNL  since  the  10-dB-down  are 
measured  from  l~^mQx  or  PNLTmax,  respectively.  Levels  in  parentheses  are 
"incomplete"  levels  in  the  sense  that  the  10-dB-down  point  was  not  reached  within 
the  measured  interval.  Appendix  D  describes  how  the  data  was  treated  before  the 
noise  metric  calculations. 

In  the  identification  section  of  Table  1-7,  the  aircraft  type  is  followed 
by  a  letter  which  uniquely  identifies  a  flight  for  which  several  simultaneous 
measurements  at  various  microphone  positions  were  made.  The  data  in  Appendix  A 
are  organized  in  the  same  sequence  as  in  Table  1-7.  The  "sentinel"  uniquely 
identifies  a  flyover  spectral  time  history  (STH)  for  computerized  retrieval  from  a 
databank.  This  sentinel,  prefixed  by  "90000,"  appears  on  the  STH  data  plots  of 
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Measured  Noise  Data:  Comments,  Waviness,  Noise  Levels 


1-73 


Identification 


Table  1-7  (Continued) 


Table  1-7  (ConHnued) 


Table  1-7  (ConHnued) 


Table  1-7  (Continued) 


designates  sound  exposure  level,  values  in  parentheses  are  derived  from  spectral  time  histories  too 
short  to  include  the  full  10-dB-down  interval. 


Appendix  A.  Some  approach  flights  occurred  displaced  laterally  by  about  230  m 
(750  ft)  from  the  overhead  position.  These  are  identified  as  "sideline"  in  Table  1-7. 
These  flights  are  not  representative  of  the  FAR  Part  36  approach  measurement 
conditions,  but  they  are  included  as  a  matter  of  interest. 

Discussion  of  Data  in  Table  1-7 

As  must  be  expected,  no  pseudotones  are  discernible  in  any  STH  for  the 
ground  (0  m)  microphone.  For  the  pole  (10  m)  microphone  pseudotone  interference 
patterns  are  sometimes  recognizable  in  the  lowermost  ! /3-octave  bands.  For 
almost  all  1.2  m  microphone  positions,  the  pseudotone  patterns  are  readily 
discernible.  From  this  first  overall  impression,  one  might  be  tempted  to  prefer  the 
ground  over  the  pole  microphone  position. 

Figure  1-27  shows  histograms  of  the  waviness  (the  mean  over  the 
10-dB-down  interval  from  Table  1-7).  Ground  and  pole  microphone  positions 
generally  show  low  waviness,  whereas  the  1.2  m  position  exhibits  generally  greater 
wavinesses  which  are  also  more  spread  out.  This  confirms  in  a  numerical  way  the 
qualitative  observations  of  the  previous  paragraph  including  the  residual  pseudo¬ 
tones  of  the  pole  microphone  position:  its  wavinesses  cluster  slightly  higher  than 
those  of  the  ground  position. 

Figure  1-28  shows  histograms  of  the  EPNL  differences  between  the 
mif  rophone  positions.  A  clue  as  to  a  preferred  choice  for  the  microphone  height  to 
minimize  pseudotones  is  not  apparent  except  perhaps  that  the  EPNL  differences 
between  the  0  m  and  1.2  m,  and  the  0  m  and  10  m  positions  do  not  cluster  around  3 
dB  as  one  might  expect  at  least  for  hard  ground,  but  around  I  or  2  dB.  This  may 
well  be  due  to  the  various  causes  of  high  frequency  anomalies  as  observed  by  Smith 
for  the  ground  microphone  (see  Section  1.4. 3. 2)  and  as  predicted  analytically  in 
Figure  1-17.  This  can  also  be  observed  in  some  of  the  STH  plots  of  Appendix  A. 
For  example,  compare  Figures  A-GI  and  A-G3:  the  0  m  data  is  above  the  10  m 
data  for  the  lower  part  of  the  spectrum,  but,  as  expected  according  to  Figure  1-17, 
crosses  over  for  the  higher  frequency  portion,  especially  close  to  overhead. 

In  summary,  while  a  complete  explanation  for  all  of  the  anomalies  at 
high  frequencies  for  the  surface  microphone  is  not  available,  it  does  appear  that 
the  basic  ground  reflection  process  resulting  from  the  Kz-inch  elevation  of  the 
ground  microphone  is  responsible  for  a  major  part  of  the  phenomena. 


Figure  1-27.  Waviness  Histograms 
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Kecommended  Correction  Procedure:  Pole  Microphone 

Table  1-8  lists  advantages  and  disadvantages  of  the  two  instrumentation 

correction  procedures  considered  in  this  study.  Most  of  those  points  have  been 

2  f 

discussed  in  previous  paragraphs.  We  were  able  to  support  Smith's  finding  that 
pole  mounting  does  not  increase  background/wind  noise.  A  cursory  analysis  of 
several  ambient  noise  records  taken  during  the  LAX  measurements  (Appendix  A) 
showed  differences  (L.grouncj  -  Lp0|e)  ranging  between  -I  and  4.9  dB  for  A- 
weighting,  and  -0.8  and  5.1  dB  for  the  overall  level.  Presumably,  this  is  true  only 
as  long  as  the  maximum  wind  speed  at  the  pole  microphone  is  below  some  threshold 
speed,  probably  the  usual  10  knots  (5  m/sec)  with  the  use  of  a  windscreen. 

We  judge  that  the  pole  microphone  position  is  preferable  to  the  ground 
microphone  position  for  aircraft  noise  certifiction  measurements.  The  only  serious 
drawback  of  the  pole  microphone  is  that  reflection  interference  patterns  (pseudo¬ 
tones)  do  appear  on  sideline  measurements.  The  ground  microphone,  however,  is 
equally  undesirable  for  sideline  measurements  as  the  sound  must  travel  close  to  the 
ground  over  much  larger  distances  than  for  the  overhead  measurements.  The  sound 
is  thus  subjected  to  both  refraction,  scattering  and  excess  ground  attenuation 
effects  for  sideline  positions  which  cannot  be  satisfactorily  removed  with  the 
present  state  of  the  art. 

We  do  not  believe  that  the  use  of  a  10  m  pole  for  microphone  mounting 
constitutes  an  undue  burden  for  an  applicant  for  an  aircraft  type  certificate,  at 
least  as  far  as  large  scale  manufacturers  are  concerned.  In  the  case  of  small 
manufacturers  or  private  hobbyists,  FAA  may  consider  making  such  poles  and 
instrumentation  available  in  some  of  the  regional  offices. 

27 

Although  Smith  suggests  o  correction  method  for  10  m  measurements 
taken  over  soft  ground,  only  hard  ground  should  be  used  for  certification  measure¬ 
ments  at  this  time.  More  research  is  required  to  ascertain  a  reliable  soft  ground 
correction  (which  depends  on  frequency)  as  a  function  of  ground  surface  imped¬ 
ance.  Such  research  would  first  have  to  start  with  methods  for  determining  ground 
surface  impedance. 
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Table  1-8 

Comparison  of  Ground  and  Pole  Microphone  Positions 


Advantages 

Disadvantages 

Ground 

No  reflection  interference 

Susceptible  to  micrometeorological 

Microphone 

whatsoever  if  flush- 
mounted. 

effects  (high  frequency  loss). 

Reflection  interference 

Flush  mounting  is  awkward  (setup 

negligible  if  mounted 
inverted  close  enough  to 

and  calibration). 

surface. 

Must  have  very  large  reflecting 
surface  to  avoid  edge  diffraction 

Easy  access  if  mounted 
inverted  close  to  surface. 

effects. 

Pole 

Reflection  interference 

On  sideline,  reflection  inter- 

Microphone 

usually  negligible. 

ference  may  not  be  negligible. 

Measurement  System  more 

Setup  and  calibration  are  cumber- 

easily  portable  than 
ground  microphone  with 
its  reflecting  surface. 

Microphone  usually  just 
above  atmospheric 
boundary  layer. 

No  increased  background/ 
wind  noise  as  compared 
to  0  m  and  1.2  m  positions. 

some. 
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1.5 


Recommended  FAR  Port  36  Actions 


Simultaneously  with  the  work  on  this  study,  Wyle  Laboratories  was  also 

working  on  other  aspects  of  aircraft  noise  certification  under  different  tasks  of  the 
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same  contract  and  under  a  different  contract.  The  recommendations  put 
forward  in  this  section  should  therefore  be  viewed  as  part  of  a  larger  set  of 
recommended  modifications  to  FAR  Part  36. 

We  recommend  that: 

o  The  current  microphone  distance  from  the  ground  of  1.2  m  be 
changed  to  10  m,  with  a  tolerance  of  plus  or  minus  0.2  m. 

o  The  800  Hz  cutoff  for  tone  correction  calculations  be  eliminated 
except  for  sideline  measurements  with  a  10  m  microphone. 

o  The  ground  surface  condition  be  more  rigorously  controlled  by 
requiring  that  the  surface  be  acoustically  hard  such  as  provided  by 
concrete  or  asphalt.  The  extent  of  such  a  surface  should  be  at 
least  25  m  (82  ft)  to  either  side  of  the  microphone  vertical 
projection  on  the  ground  in  the  flight  direction,  and  10  m  (33  ft)  to 
either  side  perpendicular  to  the  flight  direction. 

o  The  10  m  microphone  support  be  of  light  weight  slender  design  so 
as  to  minimize  its  interference  with  the  sound  field.  This  can  be 
achieved  for  example  by  using  one  vertical  pipe  (or  several  single 
vertical  pipes  stacked  end  to  end)  held  in  place  by  guy  wires. 

It  would  also  be  desirable  that  the  ground  surface  be  of  a  light  color  so 
as  to  minimize  solar  heating  effects.  However,  in  the  absence  of  more  specific 
evidence  of  the  need  for  this  refinement, it  is  not  put  forward  as  a 
recommendat  ion. 

This  change  of  microphone  height  should  not  cause  a  change  in  EPNL 

standards  (aircraft  noise  limits),  nor  should  there  be  a  correction  necessary  to 

23 

derive  EPNL  values  measured  at  1.2  m  from  those  measured  at  10  m.  Smith 
finds  that  the  two  positions  ". . .  produce  virtually  the  same  absolute  levels  over 
varied  surfaces  . . The  limited  data  of  this  study  finds  that  the  difference  EPNL 
(1.2  m)  -  EPNL  (10  m)  straddles  zero,  ranging  from  -2  to  +5  dB  (including  sideline), 
and  from  -2  to  +2  (excluding  sideline)  (see  Figure  1-28).  More  data  should  be 
examined  to  ascertain  this  point. 
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APPENDIX  A 


AIRCRAFT  NOISE  DATA  ACQUISITION  (LAX) 

Meosurement  sites  in  the  vicinity  of  Los  Angeles  International  Airport  (LAX) 
were  used  to  obtain  samples  of  measured  noise  for  different  types  of  aircraft  at 
locations  representative  of  FAR  Part  36  (takeoff,  sideline  and  approach).  The 
locations  are  shown  in  Figure  A- 1  and  described  in  Table  A-l.  The  takeoff  location 
is  much  closer  to  brake  release  than  specified  in  the  Regulation  because  of  airport 
geography  restrictions. 

Microphones  were  placed  at  ground  level,  1.2  m  and  10.0  m  as  shown  in 
Figure  A-2  over  both  a  hard  surface  (concrete/asphalt)  and  a  soft  surface 
(sand/grass)  with  the  two  surfaces  separated  by  about  10  m.  In  all  cases,  the 
microphones  were  oriented  with  the  diaphragm  horizontal.  Measurement  sites 
were  selected  to  be  in  open,  generally  flat  areas  as  free  as  possible  of  nearby 
reflecting  obstacles. 

Three  two-channel  Nagra  tape  recorders  were  used  to  record  the  acoustic 
pressure  signals  from  the  six  microphones.  Correlation  between  microphones  was 
maintained  by  using  a  common  Irig-B  time  code  generator  as  shown  in  Figure  A- 3. 
The  Irig  time  was  noted  at  aircraft  overhead  (or  sideline)  together  with 
airline/flight  number/  aircraft  type  information.  A  photograph  was  taken  at  the 
aircraft  overhead  position  (or  point  of  closest  approach  for  sideline).  The 
microphones  were  B&K  half  inch  condenser  type  4133.  Windscreens  were  used  on 
each  microphone.  Calibrations  were  performed  on  each  tape  immediately  before 
the  start  of  measurement  using  pistonphones  (B&K  types  4220  and  4230)  and  a  pink 
noise  generator  (General  Radio  1382).  The  same  three  calibrations  were  recorded 
at  the  end  of  each  tape. 

The  tapes  were  reviewed  in  the  laboratory  by  examining  the  L^time  history 
(see  Figure  A-4),  and  selections  made  for  digitizing.  Aircraft  powered  by  different 
power  plants  (high  by-pass  and  low  by-pass  engines)  were  selected  for  each  of  the 
locations. 

Digitization  of  the  selected  flyovers  was  performed  with  the  use  of  a  GR 
1926  multichannel  rms  detector  and  a  GR  1925  one-third  octave  band  filter  system. 
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Table  A-l 

Microphone  Locations  and  Surface  Conditions 


1 .  Sideline  Location  (9000  Feet  from  Brake  Release,  1600  Feet  Sideline) 

Microphone  Altitude/Surface: 

10  Meters/ Sandy  Grassland 
10  Meters/ Asphalt 

1.2  Meters/ Sandy  Grassland 

1.2  Meters/Asphalt 

*Ground/4‘  x4'x  3/4"  Plywood  Board  over  Grass 
*Ground/Aspha  1 1 

2.  Takeoff  Location  (11, 100  Feet  from  Brake  Release,  200  Feet  Sideline) 

10  10  Meters/Concrete 

1.2  Meters/Concrete 
*Ground/Concrete 

3.  Approach  Location  (7000  Feet  from  Threshold,  Under  Flight  Path) 

10  Meters/Grass  (Short  Cut) 

10  Meters/Asphalt 

1.2  Meters/Grass  (Short  Cut) 

1 .2  Meters/Asphalt 

*Ground/4'  x4'x  3/4"  Plywood  Board  over  Grass 
"Ground/Asphalt 

•Ground  microphones  were  inverted  with  1/2  inch  space  between 
diaphragm  and  ground  surface  or  wooden  board. 
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Figure  A-2.  Three  Microphone  Array  (10  m,  1.2  m  and  Ground  Microphone 
for  Measurement  of  Aircraft  Flyover  Noise  (Hard  and  Soft 
Ground  Surface). 

Note:  Ground  microphone  was  inverted,  not  embedded  into 
ground  as  shown. 
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Figure  A-3.  Instrumentation  Block  Diagram 
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Figure  A-4.  Flow  Diagram  of  a  Preliminary  Data  Reduction 
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The  average  sound  pressure  level  for  each  half  second  of  each  selected  flyover  for 
each  of  the  24  frequency  bands  was  obtained  and  stored  without  consideration  of 
time  constant  requirements,  i.e.,  no  temporal  weighted  averaging  was  incorporated 
in  the  digitization  process.  A  PDP  1 1  computer  was  used  to  write  the  digitized 
data  on  tape  and,  after  reformatting,  the  data  was  stored  by  Wyle  Laboratories  on 
a  Univac  1 108  computer  for  analysis  (one-third  octave  band  sound  pressure  level 
data  half  second  time  histories  referred  to  as  "spectral  time  histories"  or  STH). 

In  order  to  capture  where  the  aircraft  was  located  relative  to  the  microphone 
at  the  time  of  noise  emission  corresponding  to  each  half  second  of  recorded  data, 
data  from  the  FAA  ARTS  (Area  Radar  Tracking  System)  were  provided  (by 
courtesy  of  FAA).  The  aircraft  overhead  time  was  used  to  correlate  ARTS  time 
and  the  Irig-B  time  code.  The  ARTS  data  was  used  to  determine  aircraft  speed  and 
flight  path  gradient  as  well  as  distance  to  the  microphone. 

Meteorological  data  was  recorded  before  and  after  each  period  of  noise 
measurement. 

Table  1-7  in  Section  l. 4.3.3  list  the  aircraft  flyovers  selected  for  digitizing. 
The  remainder  of  this  Appendix  consists  of  spectral  time  history  plots  for  the 
selected  flyovers.  Spectra  are  shown  for  the  time  interval  approximately  between 
the  10  dB  down  points  at  I  second  intervals  (actual  data  was  obtained  at  1/2  second 
intervals).  The  sequence  of  these  plots  is  the  same  as  in  Table  1-7  (i.e.,  not 
ordered  by  dataset  sentinel),  and  is  given  again  in  Table  A-2. 

Prior  to  calculation  of  the  various  single  event  metrics,  the  data  were 
"cleaned,"  as  defined  in  Appendix  D,  by  applying  temporal  smoothing  and  correc¬ 
tions  for  ambient  noise  levels. 
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Table  A-2 


Microphone 

Location 

Approach 

Approach 

Approach 

Approach 

Approach 

Approach 

Sideline 

Sideline 

Sideline 

Sideline 


Contents  of  Remainder  of  Appendix  A 
Spectral  Time  History  Plots) 


Flight 

Sentinel 

(-9000000000) 

Page 

DC-IO/A 

15010 

A- 10 

16010 

A-l  1 

17010 

A- 12 

14010 

A- 13 

1 1010 

A- 14 

12010 

A-15 

727/B 

22010 

A- 16 

23010 

A- 17 

19010 

A- 18 

21010 

A- 19 

18010 

A-20 

20010 

A-21 

707/H 

53010 
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57010 
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Table  A-2  (Continued) 


Sideline 

B727/L 

Sideline 

B707/N 

Sideline 

DCIO/Q 

Sideline 

DC10/R 

Takeoff 

Takeoff 

727/G 

727/G 

T  akeof f 

B727/M 

T  akeof f 

B707/O 

Takeoff 

B707/P 

Takeoff 

DC- 1 0/S 

Takeoff 

DC-IU/T 

T  akeof f 

B747/U 

Takeoff 

B747/V 

68020 

A-55 

72020 

A-57 

64020 

A-58 

60020 

A-59 

69020 

A-6 1 

73020 

A-62 

65020 

A-64 

61020 

A-66 

70020 

A-67 

74020 

A-69 

66020 

A-7 1 

62020 

A-73 

71020 

A-75 

75020 

A-76 

67020 

A-77 

63020 

A-78 

34010 

A-79 

36010 

A-80 
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A-81 
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A-82 
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79020 
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A-89 
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A-93 
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A-94 

97020 

A-95 
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A- 101 

86020 
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APPENDIX  B 


SOLUTION  OF  AIRCRAFT  EMISSION  LOCATION  EQUATION 


In  Section  2.2.2. 1,  "Step  3"  requires  the  solution  of  the  following  vector 


equation: 

A  =  E  +  u  M  |  E  | 

(B-l) 

where  E  is  the  unknown  vector.  Rewrite  eq.  (B-l )  and  decompose: 

A,  -E,  =  u,  M(E,2  +  E22  +  E3¥ 

(B-2) 

A2-E2  =  u2M(E(2  +  E22  +  E32)fe 

(B-3) 

A3-E3=u3M(E,2  +  E22  +  E32)fe 

(B-4) 

Divide  (B-2)  by  (B-3)  and  (B-4),  and  isolate  E2  and  E-j,  respectively: 


E2  =  B2  +  Ej  U2/U1  ,  E3  =  B3  +  Ej  uj/uj 

where: 

B2  —  A2  *  A  |  u2/u  |  1  Bj  —  A3  —  A  |  U3/U  | 

Insert  this  into  (B-2): 

if 

A,  -E,  =U,  M  (E , 2  +  (B2  +  E ,  u2/u , )2  +  (B3  +  E ,  u3/u,)2) 


which  is  solved  for  E  j  as  given  in  Section  1 .3.2. 1 .  In  determining  which  sign  of  the 
square  root  is  appropriate  we  are  guided  by  the  fact  that  the  sound  emission  must 
occur  in  time  before  the  aircraft  is  at  the  apparent  location.  Considering  the 
extremely  simplified  case  of  A  =  (0,  0,  I),  u  =  (I,  0,  0),  then  82  =  0,  B3  =  I,  Cj  =  0, 
C2  =  -M2,  and: 


Ej  =  +  M/(l-M2)fe 

Clearly,  Ej  must  be  negative  here  to  satisfy  the  emission  precedence  conditions. 
Conversely,  had  we  chosen  u  =  (-1,  0,  0),  E|  would  have  to  be  positive.  In  general, 
the  root  has  the  opposite  sign  of  Uj. 
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APPENDIX  C 


CALCULATION  PROCEDURE  FOR  A  COMPLEX  FUNCTION 
REQUIRED  IN  THE  REFLECTION  CORRECTION 


Reference  24  shows  that  the  following  complex  function  is  required  in  the 
course  of  evaluating  the  reflection  effects  from  a  finite  impedance  surface: 

F(t)  =  I  -  TTfe  t  W(it) 


where  both  t  and  F  are  complex  numbers  and  W(z)  is  the  complex  error  function 
defined  by: 


W(z)  =  e”z  erfc(-iz) 


For  small  arguments  t  the  following  power  series  expansion  for  F(t)  is  used: 

oo 


A  x  J 


(2t2)n 


F(t)  =  I  -  TT“  t  e'  .2t‘  £  |. 

n  =  0 


For  large  arguments  t  the  following  asymptotic  expansion  for  F(t)  is  used: 

2 

F(t)  =  -2*^  U(-Re  t)  tef  +  (2t2)" 1  -  3(2t2)'2  +  I5(2t2)'3  -  ♦  . . . 

U(s)  denotes  the  unit  step  function  which  equals  0  for  negative  s,  equals  0.5  for  s  = 
0,  and  equals  I  for  positive  s. 

A  limited  numerical  study  showed  that  an  appropriate  change-over  point  from 
the  power  series  to  the  asymptotic  expansion  is  around  |  t  |  =  2. 

A  FORTRAN  function  subroutine  F0FT  which  calculates  F(t)  is  supplied  in 
Figure  C-l.  Function  values  calculated  very  close  to  |  t  j  =  2  by  the  two  methods 
described  above  are  only  approximately  equal,  not  exactly  so.  There  is  a  slight 
numerical  discontinuity  at  I  t  I  =  2. 
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COMPLEX  FUNCTION  FOFT(*,TAU> 

COMPLEX  TAU,m#T8Q, DELTA, »UM,T4 

- EQUIVALENCE  <T4,TVi) - 

IF (CAB3(TAU).LT.2. )  GOTO  914 
TSO«TAU*TAU 

- T4*WQ - 

FOFTa(Q.5,0.)/T4 

T4aT4«TSQ 

- POP  O.TS,  0  .i/T4»FQF4 - - 

T4aT«*T80 

FOFTa(1.875,0.)/T4eFOFT 

- IPCREALd-AU)  )  9 1  - 

911  DENOMal, 

GO  TO  915 

44E — PENOMaQ.S - - - 

915  FOFTaFOFT«DENOM* J« 5 44908a TAUaCEXP(TSQ) 

91)  RETURN 

*14 — TVi*W - 

TSQ*2.*TAU*TAU 

DENOMal 

- ODON*l - 

DELTAaTVC 

5UM«0 

- UERAT*0 - - 

920  SUMaSUM+DELTA 

I TERATalTERATfl 

- WW*HmU4/MX) - 

IFtTEST.LT, 5. E-5)  GOTO  910 
IF ( ITERAT .GT .25)  GOTO  9)0 

- IVE*TVE*TSa - - - 

OODNaOODN+2 

DENOMaOENOMaODDN 

- DELTAaTVE/OENQM - - 

GOTO  920 

910  TSQaTAUaTAU 

- toEW  »2.  heghsum - 

RETURN 

9)0  MR1TE(6,940)  ITERAT 

PIMMAT TIE, «  ITERATION  EXCEEDED  IN  FQFTAU  CALCULAT ION- » 1 
RETURN  1 
END 


Figure  C-l .  FORTRAN  Function  Subroutine  for  Calculating  Values  of  a  Function 
Required  in  the  Reflection  Correction. 
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APPENDIX  D 

DETAILS  OF  RAW  DATA  CLEANING  PROCEDURES 


Prior  to  calculation  of  momentary  and  single  event  metrics  (LA,  Ls,  PNL, 
PNLT,  EPNL)  the  aircraft  flyover  noise  data  measured  at  LAX  and  digitized  as 
described  in  Appendix  A,  was  subjected  to  two  processing  steps:  (I)  smoothing  each 
SPL  across  adjacent  half  second  spectra  to  provide  the  required  equipment 
response,  and  (2)  correction  for  background  ambient  noise. 

The  instantaneous  I /2-second  spectra  were  smoothed  in  the  time  domain 
according  to  FAR  Part  36,  Section  A36.3(5)  (Ref.29).  The  formula  used  in  this 
study  was: 

SPL  SPL  SPL. 

SPLSmoothed.i=IOI°9|0«)-2z  ‘  '  *  *  °'35  2  '  *  MS  z  '* 

where  SPL  is  the  sound  pressure  level  in  dB  in  any  one  I /3-octave  band,  i  indicates 
the  current  time  step  (i  -  I  is  the  previous  1/2  second,  i  -  2  is  the  previous  I  second 
before  "i"),  and  z  equals  10^’*. 

FAR  Part  36,  A36.5(d)  (3)  (Ref.29),  offers  a  choice  of  how  to  account  for 
influences  of  ambient  noise.  Either  up  to  4  bands  of  any  one  I /2-second  spectrum 
may  be  excluded  from  EPNL  calculations,  or,  SPL's  in  bands  contaminated  by 
ambient  may  be  corrected  under  an  FAA  approved  method.  Since  no  formal 
aircraft  noise  certification  was  intended  in  this  study,  it  was  decided  to  correct  the 
raw  data  (after  smoothing)  by  a  procedure  detailed  below.  We  also  recognized  that 
there  exist  two  kinds  of  ambient  noise,  one  of  acoustic  or  electrical  origin  which 
adds  to  the  signal  of  interest  on  an  energy  basis  (called  pre-detection  or 
background  noise),  and  one  of  instrumentation  origin  which  occurs  after  the  signal 
detection  in  the  spectrum  analyzer  and  does  not  add  to  the  signal  but  simply 
provides  an  analyzer  threshold  below  which  the  signal  is  masked  out  (called  post¬ 
detection  noise  floor). 

In  practice,  one  deals  with  a  mixture  of  these  two  different  noise  floor 
phenomena.  For  the  LAX  data  from  general  examination  of  many  different 
spectra,  it  was  possible  to  assume  that  there  exists  a  particular  cutoff  band  for  ail 
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spectra  below  which  the  ambient  is  energy  adding  (pre-detection)  background  noise 
and  above  which  the  ambient  is  post-detection  noise  floor.  Each  l /2-second  time- 
smoothed  spectrum  ("given")  was  compared  with  and  corrected  for  "ambient"  in  the 
following  way: 

o  Below  the  cutoff  band,  "ambient"  is  subtracted  from  "given"  on  an 
energy  basis;  however,  that  downward  correction  is  limited  to  10  dB. 

o  In  the  region  at  and  above  the  cutoff  band,  the  upper  tail  is  rolled  off  at 
a  rate  of  3  dB  per  1/3-OB  if  the  signal  gets  too  close  to  the  noise  floor. 
"Too  close"  means  that  "given"  minus  "ambient"  is  less  than  2  dB. 

A  computer  program  was  coded  which  carries  out  the  above  corrections  in 
the  following  7-step  process. 

1.  Compute  band  by  band  differences  (in  dB)  by  subtracting  the  ambient 
from  the  given.  Define  a  band  indicator  array  in  the  following  way: 

a.  All  indicators  are  first  set  to  zero; 

b.  For  those  bands  where  the  above  calculated  difference  is  less  than 
or  equal  to  a  threshold,  the  indicator  is  set  to  one.  Below  the 
pre/post-detection  cutoff  band  fN  the  threshold  is  0.457575  dB 
(this  limits  the  correction  in  step  (2)  to  10  dB).  Above  the  cutoff, 
the  threshold  is  2  dB. 

2.  For  frequency  bands  below  the  cutoff  band,  and  if  the  band  indicator  is 
0,  "ambient"  is  subtracted  from  "given"  according  to  the  energy 
subtraction  method  in  order  to  arrive  at  a  corrected  band  level: 

Corrected  Level  =  Given  +  10  log  |q  £l  -  ^"(Difference/ 1 0)j 

3.  If  either  all  indicators  are  one  (given  is  very  close  to  ambient)  or  all 
indicators  are  zero  (little  ambient  problem),  no  further  corrective 
action  is  taken  (i.e.,  skip  steps  4  through  7). 

4.  Starting  at  the  lower  end  of  the  spectrum,  a  search  is  made  for  each 
occurrence  of  a  one  indicator  for  all  bands  below  the  cutoff  band  f^. 

5.  Set  the  band  level  in  each  band  found  in  step  (4)  to  10  dB  below  the 
given  level. 
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6.  Starting  at  the  upper  end  of  the  spectrum  and  proceeding  backwards,  a 
search  is  made  for  the  first  occurrence  of  a  zero  indicator.  If  the 
highest  band  already  has  a  zero  indicator  the  next  step  (7)  is  skipped. 

7.  Set  the  level  in  the  band  just  above  the  one  found  in  step  (6)  to  3  dB 
below  the  level  of  the  latter.  Repeat  this  3  dB  per  I /3-octave  band 
rolloff  procedure  up  to  the  highest  band. 

Figures  D-l  and  D-2  show  a  spectral  time  history  at  two  steps  in  the  data 
cleaning  process:  raw  data  (O-l),  and  after  smoothing  and  correcting  for  ambient 
(D-2). 

In  most  cases,  we  were  able  to  use  the  first  spectrum  of  each  flyover  noise 
data  set  as  the  ambient  as  the  aircraft  was  far  enough  away  at  that  time.  When 
this  was  not  the  case,  an  ambient  spectrum  from  another  noise  data  set  was 
substituted  making  sure  that  microphone  location  (ground,  1.2  m,  pole)  and  data 
channel  and  approximate  time  of  day  were  the  same. 


Figure  0-1 .  Spectral  Time  History  of  Dataset  9000090020,  "Row"  Data. 
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